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Abstract 

A  theoretical  and  experimental  investigation  of  the  creep 
behavior  of  initially  flat  circular  plates  is  described  in 
which  particular  attention  is  paid  to  the  effect  on  this 
behavior  of  geometry  changes  induced  during  deformation.   A 
method  of  modeling  the  constitutive  behavior  of  a  material 
undergoing  creep  deformation  in  a  multiaxial  state  of  stress 
is  proposed.   An  energy  technique  for  solving  creep  problems 
using  this  constitutive  model  is  described  and  solutions  are 
obtained  for  the  cases  of  clamped  and  simply  supported  circular 
plates  under  a  uniform,  time  independent,  lateral  load.   The 
technique  retains  the  effects  of  both  bending  moments  and  in- 
plane  forces  in  the  deformation  of  the  plate.   A  series  of 
experiments  with  lead  plates  is  described  and  compared  with 
the  theoretical  predictions  obtained  in  this  study  and  with 
those  of  previous  investigators.   Both  the  theoretical  and 
experimental  results  show  that  the  effect  of  finite  deformation 
is  of  primary  importance  in  the  creep  behavior  of  plates  and 
must  be  taken  into  consideration  even  when  the  maximum 
deflection  is  well  below  the  plate  thickness. 
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c   =  constant  in  radial  deformation  profile  equation 

n 

e  =  extensional  strain 

h  =  plate  thickness 

K  =  material  constant  in  Norton's  law 

M  ,  M   =  radial  and  circumferential  bending  forces  per 
unit  length 

N  ,  Nn  =  radial  and  circumferential  membrane  forces  per 
r    6 

unit  length 

n  =  exponent  in  Norton's  law 

q  =  uniform  pressure  per  unit  area 

r  =  radial  coordinate  of  plate 

t  =  time 

u  =  displacement  in  direction  r  of  undeformed  plate 

V  =  internal  energy  dissipated  during  deformation 

w  =  transverse  deflection  perpendicular  to  unde- 
formed plate 

x  =  r/a 

z  =  thickness  coordinate  in  plate 

e  ,  t:„  =  radial  and  circumferential  strains 
r    0 

0  =  circumferential  coordinate  lying  in  plate 

k  ,  K_  =  radial  and  circumferential  curvatures 

r    0 

v  =  Poisson's  ratio 

a  =  stress 


w 

(-%) 

Kag 


n' 


3/4(f)2n+1 


1/n 


,  non-dimensional  parameter  used  to 
correlate  results 


(')  = 
(  )  '  = 


9t 

9_ 
3r 


(  ) 
(  ) 


INTRODUCTION 

Creep  is  a  term  commonly  applied  to  a  material  deformation 
process  in  which  the  strains  are  a  function  of  time  as  well  as 
stress.   The  phenonmenon  occurs  in  metals  when  they  are  stressed 
at  high  temperatures,  and  can  result  in  appreciable  plastic 
deformation  even  for  stresses  well  below  the  yield  stress.   As 
machinery  operating  temperatures  have  risen  steadily  over  the 
past  years  in  search  of  higher  efficiencies  the  phenomenon  of 
creep  has  been  of  increasing  interest  to  engineers.   When 
designing  components  which  will  operate  in  an  environment  where 
creep  may  occur  it  is  obviously  necessary  for  the  engineer  to 
be  able  to  predict  creep  deformations  in  order  to  produce  an 
efficient  and  workable  design.   The  mechanisms  of  creep  are 
complex  and  as  yet  imperfectly  understood,  but  enough  is  known 
to  permit  a  reasonable  analysis  of  some  cases  of  practical 
importance.   It  is  the  purpose  of  this  investigation  to  examine 
one  such  case,  that  of  a  thin  circular  plate  exposed  to  a 
uniform  lateral  pressure,  and  to  devise  a  means  of  predicting 
its  creep  behavior  including,  in  particular,  the  influence 
upon  this  behavior  of  the  geometry  changes  resulting  from  the 
deformation.   It  is  believed  that  the  general  procedure 
developed  in  the  solution  of  this  particular  problem  can  be 
applied  to  other  structural  shapes. 

Trouton  and  Rankine  observed  and  described  the  phenomenon 
of  creep  in  lead  wire  in  1904  (1).   In  1910,  Andrade  (2)  made 
the  first  systematic  investigation  of  creep  using  wires  of 


approximately  pure  lead  and  of  lead-tin  alloy.   In  1929, 
Norton  (3)  published  the  results  of  his  investigation  of  creep 
in  steel  at  elevated  temperatures  and  suggested  that  creep 
rate,  for  uniaxial  stress,  could  be  represented  as  a  function 
of  stress  raised  to  a  power  n.   This  relation,  which  has 
become  known  as  Norton's  law,  has  been  the  basis  for  much  of 
the  phenomenological  creep  analysis  in  recent  years.   Other 
functional  relations  for  uniaxial  test  data  have  been  suggested, 
such  as  a  hyperbolic  sine  function  proposed  by  Prandtl  and 
endorsed  by  Nadai  (4)  and  an  exponential  relation  suggested  by 
Soderberg  (5) .   The  greater  mathematical  tractability  of 
Norton's  law  has  made  its  use  far  more  common.   A  generalization 
of  Norton's  law  to  the  multidimensional  stress  state  was  pro- 
posed by  Odqvist  in  1933  (6) .   Odqvist  used  the  concept  of 
octahedral  stress  to  relate  the  strain  rate  components  to  the 
deviatoric  stresses.   At  about  the  same  time,  Bailey  (7) 
proposed  a  somewhat  more  complex  set  of  relations  for  general- 
izing Norton's  law  to  the  multiaxial  case  in  a  form  appropriate 
to  the  principal  directions  of  stress  and  strain.   Soderberg  (5) 
showed  how  his  exponential  creep  rule  could  be  generalized  to 
multidimensional  cases  through  the  stress  and  strain  invariants. 
Another  approach,  based  on  the  work  of  Drucker  (8)  and  Hopkins 
and  Prager  (9)  in  the  rigid  plastic  behavior  of  circular  plates, 
uses  the  Tresca  yield  surface  and  the  condition  of  normality  of 
strain  rate  vectors  to  that  surface  to  extend  a  uniaxial  creep 
law  to  a  multiaxial  problem.   This  method  was  used  by  Wahl  (10) 
for  the  analysis  of  rotating  discs  and  by  Venkatraman  and  Hodge 
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(11)  for  the  bending  of  circular  plates  subject  to  creep 
deformation.   At  the  present  time  there  seems  to  be  no  firm 
basis,  experimental  or  otherwise,  for  endorsing  or  rejecting 
any  of  the  above  theories  as  a  generally  applicable  constitutive 
relation  in  creep.   In  this  paper  another  model  for  the  general- 
ization of  Norton's  law  to  multidimensional  cases  is  proposed 

in  which  principal  strain  rates  are  assumed  to  be  a  function  of 
the  nth  powers  of  the  principal  stresses. 

The  analysis  of  the  creep  behavior  of  plates  has  thus  far 
been  confined  to  a  bending  or  infinitesimal  deflection  analysis, 
in  which  in-plane  forces  are  ignored,  or  a  membrane  analysis, 
in  which  moments  are  ignored.   Odqvist  (12)  analyzed  the  bending 
of  thin  circular  plates  with  clamped  and  simply  supported  edges 
using  his  multiaxial  relations  based  on  stress  invariants. 
Venkatraman  and  Hodge  (11)  solved  the  same  bending  problems  for 
circular  plates  using  the  Tresca  criterion  technique  mentioned 
above.   Venkatraman  and  Patel  (13) (14)  used  this  same  method 
for  the  bending  of  annular  plates.   The  Tresca  technique  was 
extended  to  the  case  of  orthotropic  bending  of  circular  plates 
by  Sandaranarayanan  (15) .   A  membrane  analysis  for  circular 
plates  undergoing  creep  deformation  has  been  done  by  Odqvist 

(12)  (15),  again  using  his  stress  invariant  generalization  of 
Norton's  law.   Onat  and  Yuksel  (17)  have  suggested  a  very 
simple  and  ingenious  technique  for  the  membrane  analysis  of 
a  thin  circular  plate  in  which  it  is  assumed  that  the  plate 
deforms  through  creep  behavior  into  a  spherical  shape.   The 
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stresses  may  then  be  assumed  to  be  everywhere  equal  and  the 
uniaxial  stress-strain  rate  rule  is  applied  directly.   There 
has  apparently  been,  as  yet,  no  mathematical  analysis  of  the 
creep  behavior  of  plates  which  retains  both  bending  moments 
and  in-plate  forces  and  which  considers  the  influence  of  the 
geometry  changes  induced  by  the  deformation.   None  of  the  plate 
analyses  cited  above  were  accompanied  by  any  experimental 
verification  nor  have  any  creep  tests  on  plates  been  noted  by 
the  author  in  the  literature.   Jones  (18) (19)  has  shown  the 
very  important  influence  of  finite  deformations  in  the  static 
and  dynamic  behavior  of  rigid-plastic  plates  and  it  would 
seem  likely  that  creep  behavior  would  exhibit  a  similar 
dependence. 

In  this  investigation  a  series  of  initially  flat  lead 
plates  of  three  different  thicknesses  is  subjected  to  various 
lateral  hydrostatic  pressures  in  a  test  apparatus  which  permits 
measurement  of  plate  deflection  at  the  center  and  at  regular 
intervals  along  the  radius.   By  observing  deflections  at 
various  time  intervals  one  is  able  to  determine  the  deflection 
rate  of  the  plate  and  the  plate's  deformation  profile.   Tensile 
creep  tests  are  also  conducted  on  specimens  made  from  the  same 
load  sheets  as  the  circular  plates  in  order  to  determine  basic 
material  creep  parameters.   Experimental  procedures,  apparatus, 
and  results  are  described  in  Chapter  I.   A  method  for  mathe- 
matically analyzing  the  creep  behavior  of  circular  plates 
including  the  effect  of  finite  deformations  is  described  in 
Chapter  II.   The  method  is  based  upon  computing  the  rate  of 
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internal  energy  dissipation  in  the  plate  and  equating  it  to  the  ex- 
ternal work  rate.  It  is  essentially  an  adaptation  of  a  technique 
suggested  by  Timoshenko  (20)  for  computing  the  large  deflection 
behavior  of  circular  elastic  plates.   The  problem  under  study 
has  several  inherent  nonlinear  aspects:   nonlinear  material 
constitutive  relations,  nonlinear  geometry  with  finite 
deflection,  and  nonlinear  change  in  creep  parameters  with  time. 
Despite  these  multiple  nonlinearities ,  by  using  the  technique 
herein  described,  the  results  for  a  circular  plate  may  be 
presented  in  the  simple  form  of  an  algebraic  equation  and  a 
set  of  graphs.   All  necessary  computer  work  is  incorporated  in 
the  graphs  from  which  the  deformation  rates  may  be  computed 
directly,  with  the  aid  of  a  slide  rule,  for  plates  made  from 
any  material  whose  creep  characteristics  are  adequately 
described  by  Norton's  law.   In  Chapter  III  the  theory  is 
compared  with  the  experimental  results  and  with  the  bending 
and  membrane  analyses  of  other  investigators. 

Both  the  theoretical  and  experimental  results  demonstrate 
the  strong  influence  of  finite  deflections  on  the  creep  behavior 
of  circular  plates.   This  influence  is  significant  even  for  a 
deflection  at  the  center  less  than  the  plate  thickness.   The 
proposed  method  of  analysis  yields  results  which  are  in  quite 
reasonable  agreement  with  the  experimental  data. 
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CHAPTER  I 

A.   PLATE  DEFORMATION  EXPERIMENTS 

1.   DESCRIPTION  OF  APPARATUS 

The  creep  behavior  of  a  series  of  circular  plates  of 

very  nearly  pure  lead  was  examined  in  the  apparatus  shown  in 

Figures  1  and  2.   The  plates  were  cut  from  commercially 

rolled  sheets  of  three  different  thicknesses  commonly  known 

as  4  lb.,  6  lb.,  and  8  lb.  plate.   These  thicknesses  were: 

Plate  Size  Thickness 
4  lb  .068" 

6  lb  .105" 

8  lb  .136" 

Local  variations  within  the  plates  from  the  above  thicknesses 

were  no  greater  than  ±.001".   All  circular  plate  specimens  and 

tensile  specimens  of  any  one  thickness  were  cut  from  the  same 

lead  sheet.   A  chemical  analysis  of  each  sheet  is  given  in 

Appendix  G.   Lead  was  chosen  for  these  experiments  because  of 

the  fact  that  it  creeps  at  room  temperature.   This  fortunate 

property  means  that  costly  furnaces,  temperature  control 

devices,  and  elaborate  instrumentation  which  would  be  necessary 

for  elevated  temperature  testing  are  not  required.   Although 

lead  is  not  a  commonly  used  structural  material  we  may  expect 

that  its  behavior  at  room  temperature  will  correspond  to  the 

behavior  of  more  common  structural  materials,  such  as  steel, 

at  higher  temperatures.   Moreover,  the  theoretical  analysis 

which  is  to  be  presented  is  valid  for  any  material  whose  creep 

behavior  is  governed  by  the  constitutive  equation  postualted. 

Thus,  if  agreement  between  theoretical  predictions 
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and  experiments  is  achieved  for  lead  plates,  then  the  analysis 
should  remain  valid  for  plates  which  are  made  from  other 
materials  whose  behavior  is  adequately  represented  by  the 
constitutive  equation. 

The  test  apparatus  consisted  of  a  short  section  of  steel 
pipe  with  a  blank  base  welded  to  one  end,  and  a  flange  and 
ring  on  the  other  end  which  was  used  to  clamp  the  test 
specimen.   The  flange  and  ring  were  machined  to  have  an 
inner  diameter  of  exactly  five  inches  and  square  edges  at  the 
inner  lip.   Centering  pegs  were  placed  in  the  flange  so  that 
the  ring  could  always  be  positioned  precisely  above  it  by 
fitting  the  pegs  into  holes  drilled  into  the  ring.   This  was 
necessary  as  the  inner  edge  of  the  ring  must  be  exactly  above 
the  edge  of  the  pressure  vessel  to  insure  proper  end  conditions 
The  flange  and  ring  had  serrations  machined  on  the  gripping 
surfaces.   These  serrations  were  relieved  in  both  surfaces  so 
as  not  to  damage  or  distort  the  edge  of  the  plate  sample.   The 
sample  plate  was  clamped  between  the  flange  and  the  ring  by 
eight  high  strength  steel  bolts  tightened  to  a  uniform  tension 
with  a  torque  wrench.   Six  dial  indicators  were  positioned 
above  the  plate  with  indicator  holders  attached  to  the  ring. 
These  dial  indicators  were  placed  so  as  to  record  the  defor- 
mations at  various  places  on  the  plate.   Positioning  of  the 
indicators  was  facilitated  by  the  use  of  extension  rods  of 
various  lengths  for  the  indicator  contact  points.   All  dial 
indicators  were  graduated  to  .0001"  and  were  checked  for 
accuracy  with  gage  blocks  before  commencing  the  experiments. 
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A.  Plate  Test  Apparatus 


B.  Plate  Test  Sample 
(after  deformation) 


C.  Arrangement  of  Dial  Indicators 


^IG.  1.2 


The  apparatus  was  pressurized  by  a  water  column.   The  water 
column  consisted  of  a  length  of  plastic  tubing  connected  to  a 
3"  diameter  plastic  cylinder  at  one  end  and,  through  a  valve, 
to  the  pressure  vessel  at  the  other  end.   The  use  of  a  water 
column  insured  that  the  hydrostatic  pressure  in  the  vessel 
could  be  both  accurately  measured  and  accurately  maintained 
over  the  long  period  of  time  required  for  a  creep  test.   The 
cylinder  on  the  upper  end  of  the  tube  was  for  the  purpose  of 
expanding  the  volume  of  the  upper  part  of  the  column  so  that 
small  changes  in  the  volume  of  the  pressure  vessel  from  plate 
deformation  could  not  make  an  appreciable  change  in  the  column 
height.   The  pressure  vessel  was  also  equipped  with  a  vent 
valve  and  a  pressure  gage.   The  pressure  gage  was  used  as  a 
check  to  insure  the  pressure  in  the  vessel  was  at  the  level  it 
should  be  and  not  blocked  in  some  way.   Test  pressure  data 
was  obtained  by  measuring  the  water  column  height  not  by  reading 
the  gage. 

2-   EXPERIMENTAL  PROCEDURE 

The  circular  plate  specimens,  having  been  cut  from  a 
single  sheet  of  commercially  rolled  lead,  were  flattened  with 
a  hydraulic  press  and  then  annealed  at.200°C  for  two  hours. 
After  cooling  to  room  temperature  a  sample  plate  would  be 
mounted  in  the  test  apparatus  as  shown  in  Figure  1.3.   The 
chamber  of  the  test  apparatus  and  the  pressurizing  column 
were  filled  with  water  before  putting  the  sample  in  place. 
The  sample  plate  was  marked  by  pencil  with  concentric  circles 
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from  the  center  at  1/2  inch  intervals,  and  then  placed 
carefully  on  the  upper  flange  of  the  test  apparatus.   The 
securing  ring  was  then  put  into  position  above  the  specimen 
and  secured  with  eight  bolts  tightened  to  a  uniform  tension 
with  a  torque  wrench.   The  dial  indicators,  having  been 
previously  attached  to  the  securing  ring,  were  then  positioned 
to  their  desired  locations  at  the  center  and  at  intervals 
along  the  radius.   The  assumption  was  made  that  deformation 
was  axisymmetric  so  that  readings  would  be  a  function  of 
radius  only.   Observation  of  the  tests  revealed  no  evidence 
that  this  was  invalid.   With  the  water  column  adjusted  to  a 
height  giving  the  desired  hydrostatic  pressure  the  test  is 
ready  to  begin.   Before  releasing  the  pressure,  initial 
readings  are  taken  on  all  dial  indicators.   These  are  the  flat 
plate,  or  t=0,  readings  to  which  all  subsequent  readings  are 
referred.   The  test  is  started  by  opening  the  valve  between 
the  pressure  column  and  the  pressure  vessel,  thus  admitting 
a  hydrostatic  pressure  to  the  back  of  the  plate.   Dial 
indicator  readings  and  room  temperature  are  subsequently 
recorded  at  various  time  intervals.   Plotting  these  readings 
at  any  given  time  gives  the  deformation  profile  at  that  time, 
and  plotting  the  readings  at  any  given  radius  over  all  times 
gives  the  deformation  history.   The  slope  of  the  deformation 
history  at  any  particular  time  is  the  deformation  rate  at  that 
time.   Of  particular  interest  is  the  largest  deformation  which 
occurs  at  the  center.   Two  tests  were  made  at  each  pressure 
using  a  different  plate  of  the  same  thickness  in  order  to  check 
the  reproducibility  of  the  results. 
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3.   RESULTS 

The  results  of  the  plate  deformation  experiments  are 
summarized  in  Figures  1.4  through  1.16.   Raw  data  is  listed 
in  Appendix  F. 

Figures  1.4  through  1.6  are  plots  of  the  deflection  of 
the  center  of  the  plate  as  a  function  of  time.   Two  ordinate 
scales  are  shown,  one  for  the  deflection  in  inches  and 
another  for  deflection  as  a  ratio  of  the  corresponding  plate 
thickness . 

Figures  1.7  through  1.9  also  show  center  deflection  as  a 
function  of  time,  but  here  all  deflections  are  plotted  relative 
to  that  at  t=l  hour  in  order  to  show  more  clearly  the  relative 
deflection  rates  for  different  pressures.   The  time  t=l  hour 
was  chosen  arbitrarily  as  a  reference  point.   Some  other  point 
could  have  served  equally  well. 

Figures  1.10  throuqh  1.16  show  representative  deflection 

w 
profiles  for  increasing  values  of  r — .   The  plots  are  of  trans- 
verse deflection  as  a  function  of  radius.   Those  shown  are  at 
time  t=l  hour.   Similar  plots  for  other  times  such  as  t=5  hours 
show  no  significant  differences.   It  can  be  seen  from  these 
plots  that  the  deflection  profile  is  approximately  bounded  by 
curves  of  the  forms 

w  =  w  (1-^-r)  2  1.1 

o    a 

and 

r2 
w  =  w  (l-^-z-)  1.2 

o    a 

where  w   is  deflection  at  the  plate  center.   For  tests  at  the 
o  ^ 

w 

lower  r—  ratios,  the  measured  profiles  correspond  quite  closely 

21 


1.25 


1.00 


35 F    -75 


50 


.25 


0 


.18 
.17 
.16 
.15 

.13 

i 

.12 


•  FIRST  TEST 
x    SECOND  TEST 


>-x 


-•-*-' 


TEST    7  8.8 
•-x * 


_•*— • x— 


8.0    PSI        _ 
• — x-x— 


I 


1  •" 

"    .10 

o 

5   .09 


fc" 


.#-•" 


■X» 


— X 


TEST    5  86 

•-x-»- 


6.1    PSI 
_•— • x— 


.08 
.07 
.06  - 
.05 


TEST    3 
• — •- 


4.1    PSI 

— •— •- 


*•  v-x--x — —x-x 
.04  Ax**  TEST    4 


-x— x x— 


f 


.03  — 
.02  — 


TEST  I  82 


.01    - 

x #-#x — • — • x-« x. 


3.0  PSI 

-x—  •* — • — • — 


1—  .00  L 


2  3  4  5 

TIME    (HOURS) 

FIG.   1.4 

Test  Results,  8  lb.  plate 


22 


si- 


1.75 


1.50 


1.25  - 


1.00 


.75 


.50 


.25 


•   FIRST  TEST 
x   SECOND  TEST 


.17 

.16 

.15 

.14 

.I3x- 

.12 


TEST    15816 


,#-*-* 


-x» 


— • — * 


-•-% — "X" 


8.0  PS  I 


~    .10 

o 
^    .09 


TEST    13814                                        6J  PSI 
_x-» — •-* • 


_v„x-*x*" 


_x_x-«-x- 


.08 

.07  I- 

.06 

.05 

.04 

.03 

.02 


TEST    12 

.x—x-x x- 


x x x x4jpS, 


r 


TEST    II 


► — •-•-•-• 


" 


TEST    9810 


2.39  PSI 
0  |  #«j«*-»-t-x-# — • » x x — •-#-• « 


oL        .oo 


J L 


2  3  4  5  6 

TIME    (HOURS) 

FIG.    1.5 

Test   Results      6    lb.    elate 


23 


3.25 
3.00 

- 

.22 
.23 

i      i      i  — i — 

•  FIRST  TEST 
x    SECOND  TEST 
TEST  24 
-  „*** x— x-x- x— x-x- 

£+■--•-•-•  •      •      • 

IT                   TEST  23 

i 

1                        1 

6.1  PSI  - 

— x — x 

2.75 

— 

.20 
.19, 

— • — • _ 

2.50 

— 

.18 
.17 
.16 
.15 

~                   TEST  21 

4.1  PSI  _ 

x x 

2.25 

f*                  TEST  22 

< 

2.00 

.14 

— 

1.75 

2 

**-* 

.13 
.12 

— 

1.50 

o 

.1  1 
.10 

—                    TEST   19 

1.84  PSI  _ 

1.25 

.03 
.08' 

.j^t-x-*— *—  * x 

^                   TEST  20 

— 

.00 

.07 
.06 

— 

.75 

.05 

— 

.50 

.04 
.03 

_ 

— 

.25 

.02 

TEST  17-18 

.82  PSI  _ 

.0,| 

.oo  L 

|£-ft • *        x—#              •"'                       "• 

0L 

1                        1                        1                        1 

1               1 

12  3  4  5 

TIME    (HOURS) 

FIG.  1.6 

Test   Results,    4    lb.    plate 


24 


t-\ 
w 

0) 

u 

(1) 
-p 
fd 

iH 
(1, 


h-      A 


CO 
9       MH 


O 


> 

•H 
•4-> 
(0 
H 

<D 


cJH  I  =  I  (g)  °M    01   13d  (Nl)  °M 


25 


CO 

-p 

i-H 
10 

a) 
u 

Q) 
-P 


H 

CO 

_•      m 
0 


o 


-p 
o 


> 

•H 
4-> 

fC 
r-i 

0) 
Pi 


dH   I  =  I  (D)  °M   .01     13d   (Nl)  °M 


26 


-><• + r  CD- 


COl 


CO 
COf 


CD 


—  ST 


-       CO 


OJ 


-4-> 

rH 

ft 


CO 


If) 
o 


o 
o 


UH  I  =  *  C5)  °M  01    ~13d  (Ml)    °M 


o> 


+  - 


UJ 


h- 


o 


o 

+J 
o 


x  I- 


CVJ 


0) 

> 

•H 
-P 

H 
0) 


27 


28 


.07 
.06 
.05- 


W, 


=  .32 


8  LB.    PLATE 

4.  !    PSI 

t  =  I  HR. 

TEST  4 

EXPERIMENT 

2    \2 


w  =  Wo(,_l_)< 


29 


30 


31 


w 
to  the  former  expression,  and  at  higher  r—  ratios  the  observed 

profile  seems  to  migrate  toward  the  latter  parabolic  form. 

This  seems  rather  reasonable  as  the  former  expression  corresponds 

to  the  profile  given  by  Timoshenko  (21)  for  the  elastic 

deformation  of  a  clamped,  hydrostatically  loaded,  circular 

plate  and  the  latter  is  the  corresponding  membrane  profile. 

The  fact  that  the  plate  behaves  in  this  way  was  of  particular 

significance  in  suggesting  the  method  of  analysis  described  in 

the  next  chapter. 

B.   TENSILE  CREEP  TESTS 

1.   DESCRIPTION  OF  APPARATUS 

Tensile  creep  tests  were  necessary  for  obtaining  the 
appropriate  creep  constants,  that  is,  K  and  n  in  Norton's  law 
which  may  be  written  as : 

e  =  Kan  1.3 

Tensile  test  specimens  were  prepared  from  the  same  lead  sheets 
from  which  the  plate  specimens  were  made.   The  tensile 
specimens  were  also  flattened  with  the  same  hydraulic  press 
and  annealed  for  two  hours  at  200°C  in  a  manner  similar  to 
the  circular  plates.   The  type  of  specimen  used  is  illustrated 
in  Figure  1.17D. 

The  tensile  testing  apparatus  is  shown  in  Figures  1.17 
A,  B,  C,  and  E.   It  consists  essentially  of  a  frame  in  which 
the  weighted  tensile  specimen  is  suspended  and  a  cathetometer , 
an  optical  measuring  device,  for  measuring  the  distance 
between  two  marks  inscribed  on  the  specimen.   There  is  also 
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a  modified  hydraulic  jack  used  to  facilitate  the  rather 
important  task  of  smoothly  loading  the  specimen  at  the 
commencement  of  a  test.   The  specimens  were  loaded  with 
calibrated  iron  weights  mounted  on  a  pan  suspended  from  the 
lower  ~^ecimen's  grip. 

The  frame  was  constructed  of  aluminum  I  beams  bolted 
together.   The  specimen  was  suspended  from  an  eye  mounted 
in  a  transverse  bar  centered  in  the  top  of  the  frame.   A 
flexible  braided  wire  with  an  eye  soldered  to  each  end  was 
used  to  suspend  the  specimen.   Each  of  these  eyes  had  its 
inner  edge  machined  to  a  knife  edge  to  insure  that  a  specimen 
would  hanq  freely.   A  slotted  bar,  as  shown  in  Figure  1.17C, 
was  used  to  prevent  any  rotation  of  the  specimen  during  the 
test.   This  was  necessary  to  insure  the  marks  were  always 
visible  by  the  cathetometer .   The  bar  was  placed  loosely  on 
its  supports  and  provided  only  a  very  gentle,  yet  sufficient, 
restraint  against  any  tendency  for  the  specimen  to  rotate. 

The  specimen  was  held  between  two  grips  as  shown  in 
Figure  1.17C.   The  upper  grip  was  made  of  aluminum  and  the 
lower  of  plastic.   The  use  of  plastic  was  to  minimize  the 
weight  of  the  lower  grip.   This  assisted  in  preventing  any 
bending  of  the  specimens  while  mounting  them  in  the  frame. 

The  initial  loading  of  the  specimens  was  accomplished 
with  the  aid  of  a  modified  hydraulic  jack.   Prior  to  the 
commencement  of  a  test  the  weights  and  weight  pan  were 
supported  by  a  small  hydraulic  jack  as  shown  in  Figure  1.17E. 
The  jack  is  adjusted  so  that  the  hook  on  the  weight  pan  is 
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barely  out  of  contact  with  the  ring  on  the  lower  grip.   When 
the  jack's  bypass  valve  is  released,  two  springs  draw  the 
piston  quickly  and  smoothly  down  from  the  weight  pan,  thus 
loading  the  specimen.   As  seen  in  Figure  1.17B  the  radial 
slots  in  the  weights  ride  on  a  verticle  rod.   This  rod 
prevents  the  weights  from  rotating  and  twisting  the  specimen 
but  does  not  inhibit  vertical  movement  of  the  weights. 

Measurement  of  tensile  strain  was  accomplished  with  the 
aid  of  the  optical  device  shown  in  Figure  1.17  A,  B,  and  C, 
known  as  a  cathetometer .   It  consists  of  two  telescopes 
mounted  so  that  they  can  be  moved  vertically  as  a  unit  or 
moved  apart  in  the  vertical  direction.   The  unit  had  a  vernier 
for  measuring  a  change  in  the  vertical  position  of  the  two 
telescopes  as  a  unit.   Each  telescope  was  equipped  with  a 
movable  hairline  whose  deflection  from  its  centered  position 
was  measured  by  a  vernier  graduated  to  .0001".   By  following 
the  separation  of  two  scribed  marks  initially  a  known  distance 
apart  on  the  specimen,  one  can  use  this  instrument  to  measure 
creep  elongation  and  have  no  physical  contact  with  the 
specimen  which  could  possibly  affect  the  normal  deformation 
process . 

2 .   EXPERIMENTAL  PROCEDURE 

The  tensile  specimens,  having  been  prepared  from  the  lead 
sheet  as  described  in  the  previous  section,  were  scribed  with 
two  thin  marks  one  inch  apart  on  the  narrow  section  of  the 
sample.   The  scribe  marks  were  made  by  clamping  a  one  inch 
gage  block  over  the  specimen  and  drawing  a  razor  blade  down 
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either  side  of  the  gage  block.   Great  care  was  necessary  in 
attaching  the  grips  to  the  specimen  in  order  to  avoid  bending 
the  soft  lead.   A  specially  prepared  block,  machined  to  the 
proper  dimensions,  was  found  useful  in  this  regard  for 
supporting  the  specimen  and  grips  while  they  were  being 
attached.   The  slotted  bar,  used  to  orevent  rotation  of  the 
sample,  was  placed  around  the  specimen  when  attaching  the 
grips.   This  assemblage  is  then  suspended  in  the  frame  by 
placing  the  hook  on  the  upper  grip  into  the  eye  at  the  lower 
end  of  the  braided  wire.   The  slotted  bar  rests  on  ledges 
placed  on  either  side  of  the  frame.   The  weight  pan  is  attached 
to  the  lower  grip  and  the  jack  is  raised  to  support  the  pan. 
The  required  weights  are  then  put  in  place  on  the  pan.   The 
telescopes  are  adjusted  so  that,  with  the  crosshairs  centered, 
the  crosshair  in  the  upper  telescope  is  aligned  with  the  upper 
scribe  mark  and  that  of  the  lower  telescope  with  the  lower 
scribe  mark.   The  test  is  begun  by  releasing  the  jack  and 
allowing  the  weight  to  be  applied  to  the  lower  grip.   As  the 
specimen  strains  the  scribe  marks  move  away  from  the  cross- 
hairs.  The  amount  of  elongation  can  be  measured  by  keeping 
the  centered  crosshair  in  the  upper  telescope  aligned  with 
the  upper  scribe  mark  and  moving  the  lower  telescope  crosshair 
to  realign  it  with  the  lower  scribe  mark.   The  amount  the 
lower  crosshair  is  moved  can  be  read  from  the  vernier.   The 
scribe  marks,  although  quite  thin  to  the  naked  eye,  are  much 
wider  than  the  hairline  when  viewed  through  the  telescope. 
The  accuracy  of  measurement  must  therefore  be  considered  as 
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about  ±.0002"  rather  than  within  .0001"  which  is  the 
graduation  of  the  vernier. 

By  measuring  the  movement  of  the  scribe  marks  over  a 
period  of  time  the  deformation  history  of  the  sample  under 
the  given  load  can  be  determined.   Performing  a  series  of 
these  tests  at  different  loadings  yields  a  family  of  curves 
showing  the  deformation  history  at  different  stress  levels. 
By  measuring  the  slope  of  each  curve  at  some  given  time  one 
can  obtain  a  strain  rate,  for  that  time,  at  each  stress  level 
tested.   Plotting  strain  rate  versus  stress  on  a  log-log 
plot  will,  in  the  case  of  a  material  whose  creep  properties 
conform  to  Norton's  law,  yeild  a  straight  line  from  which  the 
creep  constants  K  and  n  can  be  determined. 

3.   RESULTS 

Strain  versus  time  for  the  8  lb.,  6  lb.  and  4  lb.  plates 
is  plotted  in  Figures  1.18  and  1.19,  Figure  1.20,  and  in 
Figures  1.21  and  1.22  respectively.   The  slopes  of  these 
plots  are  measured  at  times  t=l  hour  and  t=5  hours.   These 
two  times  were  selected  as  t=l  hour  is  a  point  apparently 
well  within  the  primary  creep  region  where  deformation  rate 
is  changing  rapidly;  and,  by  constrast,  time  t=5  hours  has 
an  appreciably  smaller  and  nearly  constant  creep  rate  more 
characteristic  of  secondary  creep.   Creep  rate  versus  stress 
is  plotted  on  a  log- log  scale  in  Figures  1.23  through  1.25. 
These  data  indicate  appropriate  creep  constants  as  follows: 
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8  lb.  plate   n=3.7       K  =  5.56  x  10  14  @  t=l  hr . 

K  =  2.30  x  10~14  @  t=5  hrs 


6  lb.  plate   n=5.0      K  =  4.10  x  10  18  @  t=l  hr. 

K  =  1.70  x  10~18  @  t=5  hrs 


4  lb.  plate   n=4.0      K  =  8.30  x  10  15  @  t=l  hr. 

K  =  3.50  x  10"15  @  t=5  hrs 


It  is  interesting  to  note  the  marked  differences  in  the 
creep  characteristics  of  the  samples  cut  from  the  three 
different  sheets.   All  three  sheets  were  made  from  lead  of  the 
same  commercial  grade  and  one  might  have  expected  them  to  have 
exhibited  nearly  identical  creep  characteristics.   Such  was 
not  the  case,  however,  and  the  differences  are  too  great  to 
be  accounted  for  by  experimental  error  or  data  scatter. 
Gifkins  (22)  made  a  careful  study  of  the  effect  of  alloying 
on  the  creep  of  lead  using  lead-thallium  alloy.   He  found 
that,  in  very  dilute  alloys,  a  change  of  composition  of  as 
little  as  .01%  thallium  would  result  in  a  significant  change 
in  grain  size  and  creep  performance.   Creep  characteristics, 
therefore,  certainly  in  the  case  of  lead  alloys  and  probably 
also  for  other  materials,  seem  to  be  quite  sensitive  to  small 
changes  in  the  composition  of  the  material.   This  means  that 
calculations  based  on  assumed  creep  constants  may  be  greatly 
in  error  unless  the  constants  have  been  derived  from  tests  on 
exactly  the  same  material  as  that  being  considered. 
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CHAPTER  II 

Theoretical  Analysis 

A.   ASSUMPTIONS 

Consider  the  case  of  a  homogeneous,  incompressible, 

isotropic,  circular  plate  subjected  to  a  time-independent, 

lateral,  hydrostatic  pressure.   The  plate  is  assumed  to  be 

sufficiently  thin  for  the  problem  to  be  treated  as  plane 

stress  with  a   =0.   We  assume  the  creep  behavior  of  the 
z 

material  may  be  modeled  in  the  manner  discussed  in  Appendix  A, 

which  for  this  case  gives: 

•  _,  r  i        i  n  ii        I  n  t  ^    -, 

e,    =   K[    a  sgn   a      -    ?   crft      sgn   an]  2.1 

1  '    r '         ^        r        ^  '    0 '      3         0 

*Trriin  iiin  ,  ^    <i 

e0   =   K[    an         sgn   aA    -    j   o  sgn   a    J  2.2 

2  L  '  0  '    ^    0    " '  r '    '    r 

where  sgn  is  the  signum  function  defined  as 


sgn  x  = 


+1  x>0 
-1  x<0 


In  this  form  2.1  and  2.2  can  be  used  for  negative  as  well  as 
positive  stresses. 

B.   BASIC  EXPRESSIONS  FOR  MOMENTS  AND  IN-PLANE  FORCES 

Solving  2.1  and  2.2  for  a   and  a~  in  terms  of  the  strain 
^  r       0 

rates  we  find: 

,4  N 1/n  i'  ,i*  i 1/n      /'  ,i*  \  ^    ~> 

a   =  (tt— )  '    e  +i£n      sgn  (e    +?en)  2.3 

r    3K      '  r   0 '      ^     r    0 

,4  . 1/n  |*  ,i*  I 1/n      /*  ,i*  \  o  a 

0  =  (3K)     l£0  ^£r'     Sgn  (£0 +z£r^ 

The  extensional  strains  along  the  middle  surface  of  the 
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olate  will  be  denoted  as  e   and  e„    ,  ,.        .      ,.  ., 
*  r       8  and  the  curvature  of  the 

neutral  surface  as  k   and  Kft.   Adopting  the  Kirchoff  assumption 

that  sections  which  are  plane  before  deformation  remain  plane 

during  deformation,  the  total  strains,  e   and  e  ,  at  any  point 

along  the  thickness  of  the  plate  may  be  written  as: 

£   =  e   +  zk  2.5 

r    r     r 

cQ    =   eQ  +  zkg  2.6 

where  z  is  the  distance  from  the  middle  surface.   Notation  and 
sign  conventions  are  as  shown  in  Figure  2.1.   Taking  derivatives 
of  2.5  and  2.6  with  respect  to  time  gives  the  strain  rates: 

e=e+ZK  2.7a 

r    r     r 

£fi  =  efi  +  zL  2.7b 

*  , 

Substituting  into  2.3  and  2.4: 

,4  N 1/n I .   ,  i .   ,   , •  , i •  \ i 1/n  ~  „ 

°r  =  (3K}    'er    ^e0    z(Kr+iK0)|  7  2.8 

sgn[er+£eQ  +z(Kr+|K0)] 
, 4  . 1/n  I  .   ,  i  .   .    /  •  ,  l  •  \  i l/n  o  o 

ae  =  (Ik}   lee  +  ?er  +  ^V^l  2'9 

sgn[ee+ier+z  (Ke+2Kr)  ] 

From  this  point  conventional  notation  for  the  signum  function 
will  become  rather  onerous  and  we  shall  define  the  symbol 
|f(xx,  x2,  — -xn)|P 

as  meaning  the  absolute  value  of  the  function  f  is  raised  to 

the  power  p,  the  result  is  then  assigned  the  sign  which 

function  f  (x,  ,  x„,  x  )  assumes  when  evaluated  with  the 

l    Z  n 
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t   z  «0 
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r  =a 


appropriate  arguments  (x..  ,  x  , x  )  for  the  case  in 

question. 

Considering  a  homogeneous  thin  plate  we  may  compute 
in-plane  equivalent  forces  per  unit  length,  N   and  N  ,  as: 


N   = 
r 


Ne  " 


f- 

2  a  dz  2.10 

k    r 

\ 
2  a.dz  2.11 

2 


h 

M       = 
r             J 

'2 

a    zdz 

h      r 
2 

Me  =     j 

j 

rh 

2   aQzdz 
h      6 
2 

Equivalent  moments  per  unit  length  are: 

2.12 

2.13 

By  substituting  (2.8)  and  (2.9)  into  the  above  equations, 
(2.10)  through  (2.13),  and  integrating  over  z,  we  may  obtain 
expressions  for  the  plate  force  and  moment  resultants,  N  , 

NA,  M  ,  and  M„.   Some  care,  however,  must  be  exercised  in 

Or       o 

performing  this  integration.   The  signs  of  various  terms  in 
the  result  depend  upon  the  relative  magnitude  and  directional 
sense  of  the  extensional  strain  rate  and  the  curvature  rates 
in  the  plate.   It  can  be  shown,  however,  that  upon  integration 
all  cases  may  be  represented  by  expressions  of  the  form: 
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.4  .  1/n   n      1       ,i.  ,i.   ,h/.  ,i.  . i 1+1/n 
Nr  =  (3K}     H+I  .   ,  l [ler+?e9  +2(Kr+iK8)  I 

Kr  +  7Ke 


,  i •   h / .  . i •  v  1 1+1/n,  _  .  . 

er+'ee"2(Kr+iKe)l      ]  2'14 


_  ,4  .1/n   n 


N,  =  fe) 


0    X3K'     n+ 


1      r I .   , i .     h  ...  i .  x  1 1+1/n 
I  .  L  ,.    ['ee  +7er  +  2 (K6+'Kr) I 


Kfl+  ?K 

9     r 


h  , .   i .  .i 1+1/n 


I'll*  U/«,l'\lJ-"r'-'-/ilT  -.   -I  r- 


M   =  (— )1/n  -£-    X     r*l  U  +ie  +^(k  +±k  )|1+1/n 
r    (3K>     n+1  .   ,,.   l2  ] er+?e0+2 (  r  ?K0 } ' 

Kr     +iKe 


n  i.     ,i-       . h ,.     . i .    .  I 2+1/n 

e   +£eQ    +7r(«    +?Kn) 
1    r        9      2      r        0    '  s 


(2n+l)  (<r+HQ) 


HeA(£+K«)\1+1/n   + 


+   —  i  e 

2   '     r    ^    0    2V    r 


(2n+l)  (Kr+HQ) 


\K   +   *««    -   £<*+**„>  |?+1/n]  2.16 


0         2V    r 


%  (3K>  n+1    .    A    i.  l2|e0    +*WK9    *Kr'  ' 

V    *Kr 

n ..         i.  h».       i.    ,  1 2+1/n 

/*>_.t\#»,1«»  0  IT  20  IT        S 

(2n+l)  d<0+iKr) 
+  !lV^r-!'<e+^r»|1+1/n+  " 


(2n+l)  (kq+*k J 


r 


\*a+*K    "    T(<«+^)|«+1/n3  2'17 


2X"0    z"r' 's 
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C.   DEFORMATION  PROFILE  AND  STRAIN  RATES 

For  a  circular  plate  undergoing  finite  deformation  the 
curvatures  and  extensional  strains  can  be  written: 

K       =  -w"  Kfl  =  --'  2.18 

r  9     r 

e  =  „•  +  i^-2  efl  =  a  2.19 

r         2  8    r 

where  w  and  u  are  the  transverse  and  radial  deflections 
respectively  with  directional  sense  as  shown  in  Figure  2.1. 
Taking  derivatives  with  respect  to  time  we  obtain  the  rates: 

k   =  -w"  Kn  =  -—  2.20 

r  6     r 

e   =  u'  +  w'w'  eA  =  —  2.21 

r  0    r 

To  proceed  further,  the  transverse  and  radial  deformation 
rate  profiles  for  a  plate  with  a  given  set  of  edge  conditions 
must  be  determined.   For  the  finite  linear  elastic  deflections 
of  a  clamped  circular  plate  Timoshenko  and  Woinowsky-Kreiger  (23) 
employed  the  elastic  bending  analysis  profile: 

r-2 

w  =  w  (1-^-r)  2  2.22 

o    a 

where  w   is  the  deflection  at  the  center  and  a  is  the  outer 
o 

radius.   In  the  experiments  described  in  Chapter  I  it  was 

noted,  for  the  r-  ratios  examined,  that  the  measured  deformation 
h 

profile  could  be  represented  by  equation  (2.22),  especially  at 

w 
lower  values  of  r— -.   Assuming,  therefore,  that  the  vertical 

deformation  profile  may  be  represented  by  (2. 22) the  curvatures 

and  curvature  rates  may  be  written: 

4w      .,  2  4w      2 

53 


4w      -  2  4w      2 


For  the  deflection  rate  in  the  radial  direction,  u,  we 

adopt  an  expression  suggested  by  Timoshenko  and  Woinowsky- 

Kreiger  (24)  for  large  deflections  of  clamped  circular 

elastic  plates  which,  adapted  for  time  dependent  problems,  is: 

u  =  r(a-r)(c,  +  c~r  +  c,r   + c   ,r  ).  2.25 

1     2      3  n-1 

This  expression  satisfies  the  end  conditions  u=0  at  r=0  and 

r=a  and  provides  for  a  wide  variety  of  possible  profiles 

through  variation  of  the  c.  values.   The  constants,  c.#  will 
^  l  l 

be  determined  by  requiring  minimum  internal  energy  dissipation 
in  the  plate.   Taking  the  first  two  terms  of  the  series  we 

have : 

2        2       3 
u  =  c,ar  -  c,r   +  c,ar   -  c~r  2.26 

from  which: 

u'  =  (a-2r)c1  +  (2ar-3r2)c2  2.27 

£  =  (a-r)Cl  +  (ar-r2)c2  2.28 

With  these  expressions  the  extensional  strain  rates,  (2.21), 

become : 

9      16w  w  r       2  j 

e   =  (a-2r)c,  +  (2ar-3rZ)cn  +  ?-° —  (l-^-)^  2.29 

r  1  2      a        a 

2 
e   =  (a-r)c,  +  (ar-r  )c?  2.30 

Examining  the  expressions  for  e   and  e„  we  are  led  by 
dimensional  considerations  to  suspect  that  c,  and  c?  are  of 
the  form: 
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WW  WW 

I  O       O  „  _       „  I  O       O  TOT 

c,    =    c'    r-  c0    =    c'       — q—  2.31 

11a  22a 

where  c'  and  ci  are  pure  numbers.   In  order  for  each  term  in 

the  sum  to  have  the  time  dimension  each  must  contain  w  .   If 

o 

e   and  e„  are  to  be  zero  at  any  point  on  the  plate  when 

w  =0,  then  each  term  should  contain  w  .   To  be  consistent 
o  o 

with  the  last  term  in  (2.29)  the  diameter,  a,  is  selected  as 

a  length  factor.   If  one  considers  extensional  energy 

dissipation  only  and  solves  analytically  the  n=l  case  for  c. 

and  c?,  it  is  seen  that  they  are  indeed  of  the  above  form. 

Defining  x  =  — ,  we  may  now  write  the  curvature  rates  and 

extensional  strain  rates  as: 

4w 

k      =  — §•  (l-3x2)  2.32 

r     az 

4w 
Ke  =  -£?  (1-x2)  2.33 

w  w 

e   =  -^  [(l-2x)c'  +  (2x-3x2)c:  +  16x2(l-x2)2]         2.34 
r     a  1  ^ 

w  w  , 

e0  =  -fV^  [(l-x)c^  +  (x-x2)c2]  2.35 


We  now  have  expressions  for  the  curvature  rates  and  exten- 
sional strain  rates  in  terms  of  geometrical  parameters  for  a 
clamped  circular  plate.   These  expressions  can  now  be  sub- 
stituted into  2.14  through  2.17,  the  force  and  moment  resultants 

w 
For  the  higher  values  of  r—  the  experiments  showed  that  the 

deformation  profile  tends  to  the  parabolic  membrane  form: 

r2 
w  =  w   (1 — T)  2.36 

o     a 
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The  problem  may  be  solved  using  this  profile  also  and,  as 
will  be  seen,  the  solution  for  the  center  deflection  rate, 

w  ,  tends  to  converge  to  the  same  value  for  the  two  profiles 

w 

as  z —  increases.   The  same  procedure  can  be  done  with  other 
h 

assumed  profiles  as  seen  in  Appendix  D  for  the  simply 
supported  case. 

D.   ENERGY  EXPRESSIONS 

The  total  internal  energy  dissipation  in  the  plate  during 
creep  deformation  may  be  expressed  as: 

j-2tt  fa 


V 


(MrKr+M0Ke+Nrer+Neee)rdrd6  2.37 


o    o 

where  transverse  shear  energy  dissipation  is  neglected.   The 
integration,  strictly  speaking,  extends  over  the  deformed 
surface  area  of  a  plate.   However,  if  attention  is  confined  to 
moderate  transverse  deflections  and  small  strains  then  the 
deformed  surface  can  be  taken  to  equal  the  original  surface 
area.   This  procedure  is  usually  adopted  for  deriving  a 
consistent  set  of  equilibrium  equations  and  strain  relations 
when  using  the  principle  of  virtual  work.  (25)   Using  the 
previously  derived  relations,  equations  2.14  through  2.17  and 
2.3  2  through  2.35,  we  may  now  write  for  the  clamped  case: 
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„  n   _  («  ,1/"  n   2+1/n  ^  1+1/n  1  |b+a|  l+l/n 
r  r     3K     n+1         a        A  '    ' 


2+Vn  +  1  |B-A|  1+1/n 


B+A 


(2n+l)  A2  l°™>  s      T  A 


Ve  =   4)1/n  HTIh2+1/n   (^?)1+1/n   [^|E+D|1+1/n 


'E+D|2+l/n  +  l|E_D|2+l/n 


(2n+l)Dz  '    ' s        D 

+  tik+ttd*  l*-<1/ni  •  u-kM 


H  e  =  (lj)1/n  -^  h2+Vn  (!o,  1+1/n   1  [|b+a|  1+1/n 
r  r     3K      n+1  az        2A   '    ' 

w  w 

-|B-A|1+1/n]  •  [(l-2x)c'   ^  +  (2x-3x2)cJ,  ^ 

w 
+  16x2(l-x2)2  ~  ]  2.40 


Ne   -  (^-)1/n  -iL_  h2+1/n  ,V  1+1/n   1_  r|E+D|l+l/n 

Nee9  "  (3K}     n+1  h       {It)  2D  l ' b+D ' 

.   ,  w  w 

-|E-n|i+i/n]  •  [(l-x)c^  ^  +  (x-x2)c^  ^-]  2.41 


where : 

A  ::  (3-7x2) 

w  w 

B=(i-  |x)c'  ^  +  (*x  -  |x2)c^  ^  +  16x2(l-x2)2  % 

D  E  (3-5x2) 

w  w  w 

EE  (f-  2x)cJ  ^  +  (2x  -  fx2)c^  ~-     +  8x2(l-x2)2  ^ 
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With  these  expressions  2.37  can  be  written  as: 

w 

V  =   2,aM|¥)1/n  -§r  h2+1/n    (-°)1+Vn 

jK      n+l  a 


fl     (fM   +  f__   +  f__   +  f._  )xdx  2.42 


where  the  quantities  f    are  all  functions  of  x,  c|  ,  c', 

a,  -L    z. 

b 

w 

j — ,  and  n.   These  symbols,  f   ,  represent  the  terms  containing 
n  a, 

b 

x  in  2.3  8  through  2.41  which  must  be  integrated  over  the  radius. 

The  result  of  this  integration  over  x  will  be  a  function  of 

w 
c'  ,  c'  ,  : — ,  and  n  which  will  be  denoted  I.   The  equation  for 

V  mav  now  be  written  as: 


V  =  2*a2(^)1/n  Jj  h2+1/n  (!°)1+1/n  I(c.,c.,!°,n)   2.43 


The  quantities  c'  and  c'  as  previously  discussed  determine 
the  radial  deformation  rate  profile..   They  must  be  chosen  such 
that  the  rate  of  energy  dissipation  in  the  plate  is  minimized. 
Clearly,  minimizinq  V  with  respect  to  c'  and  c'  is  equivalent 

to  minimizing  I  with  respect  to  c'  and  c'.   The  minimization 

w 

of  I  must  be  done  for  the  current  value  of  r —  and  n.   The 

h 

correct  values  of  c'  and  c'  for  one  set  of  values  for  the 

w 
parameters  r—  and  n  may,  in  general,  be  different  from  those 

for  some  other  set.   One  way  of  determining  the  appropriate 

w 

values  of  c!  and  c'  might  be  to  substitute  values  for  r—  and  n 
12  h 

into  an  expression  for  I  and  take  derivatives  with  respect  to 
c'  and  c' ,  set  them  equal  to  zero,  and  solve  for  the  values 
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resulting  in  a  minimum  I.   Unfortunately,  the  expression 


{'  (fMr  +  %    +    \    +  V*d* 


is  such  that  the  integral  cannot  be  evaluated  analytically 
so  there  is  no  analytical  expression  for  I  which  can  be 
differentiated  with  respect  to  c*  and  c'.   We  must,  therefore, 
evaluate  I  numerically  and  search  for  those  values  of  c'  and 
c'  which  result  in  the  miminum  value  of  I  for  a  given  case 
(w  /h,  n) .   This  was  done,  for  the  cases  considered,  with  the 
computer  programs  whose  listings  are  qiven  in  Appendix  H. 
These  programs  were  used  with  an  IBM  1130  Computing  System 
at  Massachusetts  Institute  of  Technology.   The  computer 
output  yielded  a  series  of  smooth  curves  for  regular  variations 
of  c!  and  c'.   These  curves  showed  clearly  identifiable 
minima  with  respect  to  c'  and  c'  for  the  function  I.   The 
procedure  is  illustrated  for  a  typical  case  in  Appendix  E. 

E.   SOLUTION  FOR  DEFORMATION  RATE 


If  a  structure  is  assumed  to  deform  in  a  manner  prescribed 
by  a  kinematically  admissible  virtual  velocity  field,  (i.e., 
continuous  and  satisfies  the  velocity  boundary  conditions) , 
then  the  principle  of  virtual  velocities  states  that  the  total 
virtual  internal  energy  dissipation  is  equal  to  the  total 
virtual  external  energy  dissipation. 
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w 
Having  minimized  the  integral  I  for  some  case  (r— - ,  n)  we 

find  the  deformation  rate,  w  ,  by  using  the  principle  of 

virtual  velocities  and  considering  a  small  variation  of  the 

deformation  rate  about  its  equlibrium  value.   This  gives: 

dV  ,_._  „  ra 

o 

where  q  is  a  time-independent,  transverse  pressure.   For  the 

case  of  the  clamped  plate  we  have: 

2tt  /aq6w  rdr  =  27rq6w  fa  (1-^r)  2rdr  =  ^-  6w  2.45 

o  o    o  a.  5  o 


-TTT   *5w    =   2tt  /  q<5w  rdr  2.44 


and: 

dV  x .  .  1/n  ro   2/4  xl/n,l  . 1+1/n,  2+l/nx1 x .        0  „, 

-r-r   6w  =  w   '   [27raz(-^j7)  '  (—2-)      h   '  I]6v;       2.46 

dw   o  o          3K     a                 o 
o 

The  expression  2.46  is  the  derivative  of  V  for  some  specific 

w 
value  of  r—  and  the  derivative  is  properly  taken  considering 

w 

r—  as  a  constant.   We  are  seeking  the  value  of  w   for  a  given 

w 

instant  when  r-^  has  some  fixed  value.   An  instant  later  w 
h  o 

will,  of  course,  have  some  other  value  and  at  that  time  w 

o 

will   also   have    some   other  value,    but   at   a   given    instant  w  /h 

may   be    considered    fixed.       Equating    2.45    and    2.46   we   obtain 

the    result: 

2  +  2/n 

w      =    [     ,    q?, oTw Jn  2'47 

(f-)1/n   h2+1/n    61 
3k 


or 


*o   ■    [<R>2+Vn    <^>1/n   q(|iUn    -  2.48 
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If  we  define  the  quantity  $  as: 
w 

*=   t(^}2n+lJ1/n 
3/4(a)2n+l 


This  can  be  written  in  the  non-dimensional  form: 

w  ^  ,  i 

/  o   x    .,  , .  ,a.  zn+1   .n  _  _A 

(j^n)  =  3/4  (E)       *  .  2.50 

The  parameter  $  is  a  direct  function  of  I,  the  form  of  this 
function  depending  upon  the  assumed  velocity  profile.   For  the 
clamped  case: 

4>  =  (|y>  2.51 

It  is  convenient  to  define  $  as  it  permits  a  direct  correlation 

of  results  among  the  different  velocity  profiles  considered. 

When  the  quantitv  I  has  been  minimized  with  respect  to  ci  and 

w 

c',  4>  is  a  function  of  n  and  ~.   Values  of  $ ,  for  the  clamped 
2  h  ^ 

w 

case,  are  plotted  versus  r— ■  for  various  values  of  n  in 


Figure  2.2.   Use  of  these  values  with  any  of  the  equations  2.47 
through  2.50  will  give  the  center  deflection  rate  w   as  a 
function  of  the  applied  pressure  and  the  plate's  material  and 

qeometrical  oroperties.   Once  w   has  been  determined,  the 

c  o 

deflection  rate  at  any  other  point  on  the  plate  can  be  found 
from  the  assumed  velocity  profile,  which  in  this  case  is  given 
by  equation  2.22. 
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It  was  noted  in  the  plate  experiments  that  for  the 
w 
higher  values  of  r-^-  the  profile  began  to  tend  toward  a 

parabolic  form  (or  membrane  profile)  described  by  the 

equation  2.36.   The  entire  analysis  just  described  can  be 

repeated  using  the  profile  2.36  vice  2.22  to  obtain  values 

of  $  for  the  parabolic  configuration.   When  this  is  done  it 

is  found,  as  shown  in  Figure  2.3,  that  the  difference  between 

w 
the  two  solutions  is  small  for  r—  =  1.5  and  has  virtually 

w 
disappeared  when  r-^-  =3.0.   It  seems,  therefore,  that  when 

w 
the  plate  deforms  to  a  value  of  r—   large  enough  for  there  to 

be  a  significant  deviation  from  the  profile  of  equation  2.22 

it  will  have  deformed  enough  so  that  the  solution  for  w   is 

3  o 

relatively  insensitive  to  the  assumed  profile. 

The  computer  program  listed  in  Appendix  H  is  written  for 
an  assumed  velocity  profile  of  the  form: 

w  =  w  fa.  (l-^-)2  +  a.(l-f4)l  '  2.52 

ola       2    a 

where  a,  +  a~  =  1.   The  analysis  utilizing  this  profile  is 
outlined  in  Appendix  C.   If  we  assign  the  values  a,  =  1  and 
ot_  =  0,  the  analysis  with  2.52  is  the  same  as  for  the  clamped 
prof ile, 2 . 22 .   If  we  assign  values  a,  =  0  and  a~   =  1  the 
analysis  is  for  the  parabolic  membrane  profile,  2.36.   Some 
experimentation  with  fractional  values  for  a,  and  a~  was 
attempted,  but  these  produced  no  more  useful  results.   The  best 
aqreement  with  experimental  data  came  with  a,  and  a~  at  their 
extreme  values  as  above.   The  final  result  using  this  profile 
is : 
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FIG.    2.3 

Co^oarison  of  clamoed  and  membrane  solutions 
for  n  =  4.0 
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w  o  , i   3-a, 

(-2-)  =  3/4(^)2n+1  (— L-) n  2.53 

Kaq  C 

so  in  this  case 

3-a, 
*   =  (u-H  2.54 

c 

where  I   is  the  appropriate  integral  for  the  combined  profile. 

F.   SIMPLY  SUPPORTED  CASE 

The  simply  supported  circular  plate  may  be  treated  in 
the  same  manner  as  above  if  the  velocity  profile  appropriate 
to  it  is  known.   In  the  elastic  bending  analysis  of  this  case 
Timoshenko  and  Woinowsky-Krieger  [26]  determined  a  deflection 
profile  of  the  form: 

w  =  w  (1-^r)  (1-ix^  ^t)  2.55 

o    az     5+v  az 

which  for  v  =  1/2  gives  the  velocity  profile: 

w  =  w  (i-£y)  (l-l  ^  .  2.56 

o    a      11  a 

There  is  no  experimental  data  available  in  order  to  verify 
this  profile  as  in  the  clamped  case,  but  at  least  a  reasonable 
estimate  of  the  creep  behavior  of  the  simply  supported  case 
can  be  expected  from  an  analysis  with  the  profile  of  2.56. 
This  analysis  is  outlined  in  Appendix  D  and  yields  a  final 
result  of: 

w 

(     )  =  3/4(-)      (ir_)  2.57 

Kaq  s 


65 


so    that   in   this   case: 

*  =  '22V'  2-58 

s 


where  I   is  the  integral  appropriate  to  the  simple  support 

w 

profile.   Curves  of  $  versus  . —  for  various  values  of  n  are 

h 

given  in  Figure  2.4.   In  Figure  2.5  the  results  of  the 
clamped  and  simply  supported  cases  are  compared  for  the 
case  n  =  4 . 
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CHAPTER  III 

Comparison  of  Theoretical  and  Experimental  Results  and  Results 

of  Previous  Investigations 

A  comparison  of  the  experimental  and  theoretical  results 
is  shown  in  Figures  3.1,  3.2,  and  3.3.   Here  the  non-dimensional 

w  9+1 

quantity  ( J/^  "    ^s  pl°tte<^  versus  the  non-dimensionalized 

Kaq       w 
center  deflection,  r— ,  for  the  creep  parameter  n  appropriate 

to  each  of  the  three  types  of  plate  tested.   In  these  plots 

the  appropriate  analytical  results  for  comparison  with  the 

experiments  are,  of  course,  those  for  the  clamped  case.   Also 

shown  are  plots  of  the  bending  analysis  results  of  Venkatraman 

and  Hodge  [11] ,  and  of  Odqvist  [12] ;  and  the  membrane  analyses 

of  Odqvist  [12]  [16],  and  of  Onat  and  Yiiksel  [17].   The 

mathematical  details  of  reducing  these  various  analyses  to  a 

form  appropriate  to  this  plot  are  given  in  Appendix  B.   The 

bending  analysis  solutions  are  horizontal  lines  since  that 

type  of  analysis  does  not  consider  geometry  changes.   The 

membrane  analyses  consider  geometry  changes  but  not  bending 

moments  within  the  plate.   Experimental  points  are  given  at 

each  pressure  tested  for  times  of  1  hour  and  5  hours.   These 

points  are  found  by  obtaining  the  deflection  rate  at  the 

indicated  times  from  the  plate  test  experimental  data  and 

substituting  this  value  along  with  the  appropriate  creep 

constants  found  in  the  tensile  tests  and  the  pressure  and 

w        2   , 

geometrical  constants  into  the  form  ( J/^     '      The 

Kaq 
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FIG.  3.1 

Comparison  of  experiment  and  theory, 
8  lb.  plate 


70 


c 
o  CT 

>  ° 


o| 


Comparison  of  experiment  and  theory , 
6  lb.  plate 
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Comparison  of  experiment  and  theory, 
4  lb.  plate 
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value  thus  obtained  is  then  plotted  versus  the  appropriate 

w 
value  of  r—  as  determined  from  the  plate  test  data.   The 

times  t=l  hour  and  t=5  hours  were  chosen  as  representing 

cases  of  rather  rapid  and  relatively  slower  deformation 

respectively. 

The  marked  influence  of  finite  deformation  on  the  creep 
behavior  of  a  clamped  circular  plate  is  clearly  indicated  by 
both  the  theoretical  prediction  of  equation  2.50  and  the  ex- 
perimental results,  which  are  seen  to  be  in  substantial 
agreement. 

A  similar  plot,  Figure  3.4,  shows  the  relationship  of 
the  simply  supported  case  analysis  to  the  results  of  the 
simple  support  bending  analysis  of  Venkatraman  and  Hodge  [11] 
and  of  Odqvist  [12]  and  the  membrane  solutions  of  Odqvist  [12] 
[16]  and  of  Onat  and  Yii  ksel  [17]  ,  for  the  case  n  =  4.   No 
experimental  data  is  available  for  this  type  of  end  condition, 
but  the  analytical  results  also  indicate  a  marked  influence 
from  finite  deformation. 
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FIG.  3.4 
Comparison  of  simply  supported  ^late  theories 
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Summary  and  Conclusions 

Experimental  and  theoretical  results  have  been  presented 
which  show  that  the  effect  of  finite  deformation  upon  the  creep 
behavior  of  circular  plates  is  of  primary  importance.   It  is 
necessary  to  consider  this  effect  even  when  the  maximum 
deflection  of  the  plate  might  be  well  below  the  plate  thickness. 

A  method  of  modeling  the  constitutive  behavior  of  a 
material  undergoing  creep  deformation  under  a  multiaxial  state 
of  stress  has  been  proposed.   A  technique  for  solving  multi- 
axial  creep  problems  by  using  an  energy  method  and  the 
suggested  constitutive  equations  yields  a  result  in  the  form 
of  an  equation  and  graphs  which  can  be  quite  easily  applied. 
Solutions  were  obtained  for  a  circular  plate  under  a  time 
independent  uniform  lateral  load  with  clamped  and  simply 
supported  edge  conditions.   A  comparison  of  theoretical  and 
experimental  results  for  the  clamped  case  showed  substantial 
agreement;  indicating  that  the  analytical  technique  described 
may  be  expected  to  provide  quite  reasonable  estimates  of  plate 
deformation  rate  when  the  creep  behavior  of  the  material  may 
be  adequately  represented  by  Norton's  law,  and  the  appropriate 
values  of  the  material  creep  parameters  are  known.   A  comparison 
of  these  results  with  those  of  several  previous  investigators 
has  also  been  presented.   Lead  was  selected  for  the  experimental 
material  in  these  tests  but  the  analytical  results  should 
remain  valid  for  any  material  whose  behavior  is  adequately 
represented  by  the  constitutive  equations. 

75 


Specific  consideration  was  given  only  to  circular  plates; 
it  can  be  expected,  however,  that  the  conclusion  as  to  the 
importance  of  finite  deformation  will  apply  as  well  to  plates 
of  other  shapes,  initially  straight  beams,  and  to  initially 
straight  cylindrical  shells. 
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APPENDIX  A 
Constitutive  Equations  for  Multiaxial  Creep 

Creep  data  is  not  plentiful,  and  that  which  is  available 
is  almost  always  the  result  of  uniaxial  tests.   It  has  been 
found  that  much  of  this  experimental  creep  data  can  be  repre- 
sented by  an  expression  of  the  form: 

e  =  K  an  (A.l) 

where  K  and  n  are  material  parameters  which  are,  in  general, 

functions  of  time  and  temperature.   This  expression  is 

commonly  known  as  Norton's  law.  [3] 

In  generalizing  Norton's  law  to  a  form  suitable  for 

multidimensional  problems  the  most  common  approach  has  been 

to  use  the  concept  of  equivalent  stress  and  equivalent  strain 

rate  for  determining  the  proportionality  between  strain  rate 

and  deviatoric  stress.   Mendelson  [27]  uses  a  relation  of  the 

type: 

• 

3   e 

e .   —  —  —  s 
2.J         2    a     i  j  (A.  2) 

where  creep  equivalent  strain  rate,  £  ,  may  be  related  to  the 

e 

equivalent  stress,  a  ,  by  Norton's  law. 

e 

Mere : 

o   =  /3JT  (A. 3) 

e       2 

and  e   is  the  strain  rate  resulting  from  a  stress  of  this 

e  3 

magnitude.   Odqvist  [6]  uses  the  expression: 

Eij    2  ec    (o    '  a  (A. 4) 

c       c 

• 

where  e   is  some  standard  value  of  creep  rate  for  which  a=o    . 

c  r  c 
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These  expressions  seem  reasonable  and  satisfy  the  mathematical 

requirements  for  multiaxial  creep  constitutive  relations,  but 

are  often  rather  intractible  for  analytical  work  since  one 

must  always  compute  a   in  order  to  know  the  proportionality 

between  the  local  stress  and  creep  rate. 

In  this  analysis  the  constitutive  behavior  of  a  material 

undergoing  creep  deformation  in  a  multiaxial  state  of  stress 

will  be  modeled  in  the  following  way: 

*T.rrn        n       n,  , ,  r-  \ 

c  =  K[a.   -  v  a.   -v  a,  ]  (A. 5) 

1      1       j      k 

where  the  subscripts  i,  j  ,  and  k  represent  principal  directions 
of  stress  and  strain  rate.   It  can  be  shown  that  for  conser- 
vation of  volume  we  must  have  v  =  1/2.   In  real  materials  such 
things  as  voids  and  inclusions  may  give  rise  to  situations  in 
which  volume  is  not  conserved.   If  this  effect  is  important 
then  in  such  cases  some  other  value  for  v   must  be  assigned  as 
appropriate  to  the  specific  material  under  consideration. 

It  is  also  desirable  that  the  constitutive  relation  (A. 5) 
hold  for  both  positive  and  negative  a  for  all  values  of  n. 
The  significance  of  positive  and  negative  signs  here  is  merely 
one  of  direction,  tensile  or  compressive,  and  this  directional 
sense  must  be  retained  for  all  values  of  n.   This  may  be 
accomplished  through  the  use  of  the  signum  function  defined  as: 

s^  *  -  <  -i    x<°  (A-6) 

An  expression  of  the  form  |x|   sgn  x  is  used  to  denote  that  the 
absolute  value  of  x  is  raised  to  the  n  power  and  the  result  is 
assigned  the  sign  of  x.  i 
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For  homogeneous  isotropic  materials,  therefore,  the 
proposed  constitutive  relation  takes  the  form: 

ei=K[|ai|n  sgn  a±-^  |  a^  |n  sgn  a_.-^|ak|n  sgn  aR]        (A. 7) 

These  equations  are  valid  for  a  material  which  has  the 
following  properties: 

(1)  Homogeneous,  incompressible,  isotropic  material. 

(2)  Coaxial  stress  and  strain  rate  tensors. 

(3)  Zero  deformation  rate  under  hydrostatic  pressure. 

(4)  Coordinates  correspond  to  principal  directions  of 
stress  and  strain. 

(5)  Creep  behavior  corresponds  to  Norton's  law  in 
uniaxial  stress. 

It  should  be  noted  that  this  model  predicts  an  effect  from 
a  hydrostatic  pressure  superimposed  on  a  deviatoric  stress. 
The  validity  of  this  may  be  questionable  and  was  not  really 
tested  in  the  cases  studied  in  this  investigation.   For 
applications  where  larger  hydrostatic  pressures  are  involved  it 
may  be  better  to  use  only  the  deviatoric  components  of  the 
stresses . 

For  the  case  of  a  thin  circular  plate,  choosing  cylindrical 
coordinates  and  assuming  o^    =    o      =0,  equations  A. 7  become: 

er  =  KUar|n  s9n  ar  "2|a0|n  sgn  aQ] 

e   =  K[|aQ|n  sgn  aQ  ~i|ar|   sgn  a rl  (A. 8) 

£z  =  K[-£|ar]n  sgn  ar  -i|aQ|n  sgn  aQ]. 
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APPENDIX  B 
Interpretation  of  Results  of  Previous  Investigators 

Introduction 

The  following  theoretical  treatments  all  apply  to  circular 
plates  loaded  with  a  time  independent  pressure  uniformly 
distributed  over  the  entire  plate  area.   The  infinitesimal 
solutions  in  Sections  I  and  II  retain  only  bending  moments 
while  solutions  in  III  and  IV  consider  geometry  changes  but 
retain  only  membrane  forces. 
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I .   Venkatraman  and  Hodge  Solution  [11] 
a.   Clamped  Edges 

The  maximum  deflection  rate,  at  the  plate  center,  is 
expressed  in  the  form 


•  • 


9      2 

w   =  w*  GaZ    (3^-)n  (B.I.I) 

o  12y 

where  G  is  an  arbitrary  function  of  time  and 

~,    n    ,h.  2+1/n  /Ti    T  _. 

U=  2X    2nTT   (2)      '  (B.I. 2) 

The  term  A  is  obtained  from  the  assumed  form  of  Norton's  law 
given  as 

e  =  (j)n.  (B.I. 3) 

When  compared  with  the  form  given  by  equation  A.l  which  was: 

e  =  Ko11  (B.I.  4) 

we  have 

\   =   4  Vn  (B.I. 5) 

Substituting  B.I. 2,  B.I. 3,  and  B.I. 5  into  B.I.I: 

/,*  ^,a,2n+l  „   n  ,  22+1/n (2n+l) ,n  ,D  T  n 

wq  =  w*  G(^)      Kaq   [ j—' H  (B.I.  6) 

or 

• 

w 
( 2_) 

K^ m      ^         22+Yn(2n^l)]n  ^ 

■*  i„.i  ^4n 

.a. 2n+l 

V 

which  interprets  B.I.I  in  terms  of  parameters  introduced  in 
Chauter  III. 
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For  comparison  with  the  results  in  this  paper  we  take 
G=l.   The  term  w*  is  plotted  as  a  function  of  n  in  reference 
[11]. 

b.   Simple  Support 

For  simply  supported  edges  the  deflection  rate  at  the 

plate  center  is  of  the  form: 

9     2 
w   =  w*  Ga  (-3^ — )n  (B.I. 8) 

o  6y 

Substituting  and  rearranging  as  above  this  becomes: 

w 

Kaq 


(-  -? ) 


,,2+1/n ,_  ,.  . 
2    '  (2n+l) , n 


•a. 2n+l 
V 


=  w*  G[— 19r/     ']  (B.I. 9) 

lzn 


where  as  before  w*  is  plotted  in  reference  [11] . 
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II.   Odqvist  Bending  Solution  [12] 
This  result  is  of  the  form 

w   =  a  n  n  V(o)  (B.II.l) 

o 

where  ^ 

=  2n+l   3a   2n    ga  (B.II.2) 

n     ,2      ha 
h         c 

and  V(o)  is  plotted  as  a  function  of  n  in  reference  [12]  for 

the  clamped  and  simply  supported  cases.   The  term  o   is  from 

the  assumed  constitutive  relation: 

eij    2[o    '  a  (B.II.3) 

J      c       c 

which  for  uniaxial  stress  reduces  to 

e  =  i  n   an  I  (B.II.4) 

c 

and  when  compared  with  our  notation: 

e  =  Ka11  (B.II.5) 

gives : 

°c   =  (l)Vn  (B.II.6) 

Substituting  these  expressions  into  B.II.l  and  rearranging  we 
have: 

w 

( — ) 

v      n  o  j.1        n+1 
Kag =    (in±l)n   (3)~2~~  V(o)  (B.II.7) 


,a. 2n+l 


n 
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III.   Odqyist  Membrane  Solution  [16] 

In  this  analysis  the  transverse  deflection  at  the  center 

of  the  plate  is  expressed  in  the  form: 

2a   h  W(o) 
c 

q 


w   =  9 (B.III.l) 


where 

W(o)  =   ^-    {l^C+[(4n^(^}+  9)]  C2},         (B.III.2) 


P.2  .      fA~2. 


p   -Htj.  ,  (B.III.3) 


c 


2 


£ £j -^   =  -C,  (B.III.4) 

8(n+3)s 


s  =  [ L_ —  +  8(n+2)(n+3)  — ]  n+2  f  (B.III.5) 

1\  +  Z  i 

s  p 

h  is  half  the  plate  thickness,  and  C  is  tabulated  as  a  function 
of  n  in  the  reference.   In  this  form  equation  B.III.l  gives  the 
center  deflection  as  a  function  of  time.   For  the  deflection 
rate  we  take  the  derivative  with  respect  to  time: 

2a  h 

w   =  — —    V7(o)  .  (B.III.6) 

O     q 

Differentiating  B.III.2  with  respect  to  time: 

2 

A  i    ^               Pi         s         ;,^    i_          (4n2    +    27n   +    9)C2-,  ,D    TTT    71 

W(o)    =    -    4-        -2       {1+    ic    +    i8(n+3) }  (B.HI.7) 

from    B.III.4: 

s    =    -C    8(n+3)sn+3  (B.III.8) 

P    2 
1 
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Substituting   B.III.5    and   B.III.8    into   B.III.7: 

,'w^\         ->m    xn  n^l^j.    <4n*    +27n   +   9)C2t 
W(o)    =    2C(n+3)  [1+*C+   i8(n+3) 

1  ct    -— 

[-—rro   +   8(n+2)(n+3)    —  ]    n+2  (B.III.9) 

~   n+z  2 


Substituting   B.III.5    into   B.III.2: 

p/         1  Ct      -L. 

W(o)    =   ~[      ^+2   +    8(n+2)(n+3)    ^^]n+2 


P    2 


ti^c+  (4n21;^;3;9)  c«]  (b.hi.i<» 

Solving    for   t: 

0     2 

t   =                 * r/    ^    w(o)     xn+2  _    _JL_i                  /»    ttt    11) 

fc         8C(n+2)(n+3)     l(         Tb^)  ^TT2J                  (B. III. 11) 

p  s 

1  o 


where  B   =    1   +   *C   +    (4n'    +2?"   +   9)      C2         .  (B. III. 12) 

18(n+3) 


Substituting  for  t  in  B.III.9: 

ni,  n      2n+2 

In+2  fn+3)       pi 


•w(o)  =  2C  °     (n+3) 
,n+l 

2n+2 


(w(o))n+1 


P 
i 


(W(o))n+1  (B.III.13) 

where  _ 

i  _  2C  Bn  Z    (n+3) 

4n+l  .  B. III. 14) 

Note  that  4>  here  is  that  of  reference  [16]  and  bears  no 
relation  to  the  quantity  $  as  used  in  Chapter  II. 
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G     = 

(J_,n     0» 

c 

giving: 

a 
c 

Odqvist's  constitutive  relation  was: 

J  c       c 

which  reduces  in  uniaxial  stress  to: 


(B. III. 16) 


(B. III. 17) 

Substituting  B.III.l,  B.III.3,  B.III.6,  B. III. 14,  and  B. III. 17 
into  B. III. 13  and  rearranging: 
w 

t    °     ) 

S23 =    i  (2h)n+l  (B. III. 18) 

w 

,a_.  2n+l  ° 

4h; 

Odqvist  uses  the  symbol  h  to  denote  half  of  the  plate  thickness 
rather  than  the  full  thickness.   If  we  call  the  full  thickness 
h,  we  have  a  form  analagous  to  that  of  the  other  results: 

w 

( 9~) 

__Kag_n 


=      $      ± (B. III. 19) 


h  (jj-.) 

where 

i      =      2C[l+jC  +   4n^„ffi  +   9      CM"+2    (n+3) 

.n+1 
4 

and  C  is  tabulated  as  a  function  of  n  in  reference  [16] . 
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IV.   Membrane  Solution  for  a  Circular  Plate ,  Onat  and  Yiiksel  [17] 

Onat  and  Yiiksel  assume  that  an  initially  flat  circular 
plate  deforms  into  a  spherical  surface  upon  application  of  a 
lateral  hydrostatic  pressure.   Consideration  of  this  geometry 
shows  that  the  middle  surface  strain  rate  tensors  have  prin- 
cipal directions  which  are  directed  along  the  sphere's 
parallel,  meridian,  and  the  normal  to  the  surface.   The  first 
two  strain  rates  are  positive  and  the  third  is  given  by 

e^  =  7-  =  -5—  =  -  atan   ■=•  (B.IV.l) 

3    h     da  2 

For   a   spherical   membrane: 

a      =   a «   =   a   =   -?-£-  =   -=~-     -. 2-rrr r  (B.IV.2) 

12  2h  2h        sma       (1+cosa) 

o 

Creep  deformations  are  governed  by  the  relation: 

a  tan  |  =  f (a)  (B.IV.3) 

where  for  Norton's  law 

f  (a)  =   e  (— )n  -  -^   an  =  Kan  (B.IV.4) 

o  a        n 
o     a 


The  symbol  p  in  B.IV.2  above  represents  the  uniform  pressure 

per  unit  area  as  used  in  [17] ,  D   is  the  diameter  of  the  built- 

o 

in  support,  and  h   is  the  original  membrane  thickness. 
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Figure  B.IV.l 


From  the  geometry  of  the  above  figure  we  see 
w   =  R(l  -  cos  a) 


(B.IV.5) 


=  R  sin  a 


Equating  R 


sin  a 


D 


(1-cosa) 


2  w 


(B.IV.6) 
(B.IV.7) 


.                         ot   i—cosci 
By  trigonometric  identity:   tan  -x-  =  i , 


2w 
therefore,    tan  ~   =— — . 

o 


(B.IV.8) 


Differentiating  with  respect  to  time 

,  •  2w 

dt  (tan2)  =  2   SGC  2  =  D~ 

o 


4w 


or   a   = 


D  sec^ 
o     2 


(B.IV.9) 
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Substituting  B.IV.2,  B.IV.4,  B.IV.8,  and  B.IV.9  into  B.IV.3 

• 

4w       2w  Dp  , 

/____0__\  (  °  \  -        v      I         °     )U      / i \U 

*      2aM  D    '         l   2h   ;   l  sina  (1+cosa)  ; 
D  sec  r    o  o 

o    2 

and  rearranging 

w 

( — S — > 

n 
Kap  -, 


=  [•?(— i 7T-; r)   sec*-=r] 

.  -,  w  sina  (1+cosa)        2 

(-)n+1  (— ) 

ho  ho  (B.IV.10) 

u  D 

where    a  =   o  . 


The  non-dimensional  parameter  resulting  from  the  energy 

analysis  is 

w 

( 2__) 

Kapn 


,a_v  2n+l 

V 
o 


Transforming  B.IV.10  to  this  form: 


w 

,affi+l "    .a.,/  w  ^'sina'  (1+cosa) '"    sec'f]  (B.IV.ll) 

o 

w 
Given  : —  and  — ,  angle  a  can  be  found  from  B.IV.7  and  values  of 
o      o 

• 

w 
o 


n  w 

KaP  j.j     j=        x.-      ^o    ..a 

— % — -s-  computed  as  functions  of  r —  and  r-. 

a^  2n+l  h       h 

h' 
o 


(")  o      o 
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APPENDIX  C 
Mathematical  Details  for  Analysis  with  Combined  Bending  and 
Membrane  Profile  for  the  Clamped  Circular  Plate 

Experiments  indicate  that  a  clamped  circular  plate 
undergoing  creep  deformation  from  a  transverse  hydrostatic 
pressure  assumes  a  velocity  profile  described  by: 

w  =  w  U---t)2  (C.l) 

W  W 

when  r—  is  relatively  small.   At  larqer  values  of  -. —  the 
h  ^  h 

profile  tends  to  the  membrane  form  given  by; 

w  =  w  (l-f-r)  (C.2) 

o    a 

These  two  forms  can  be  combined  by  writing  the  velocity 
profile  as: 

w  =  ^[c^U-!^-)2  +   a2(1~rz")]  (C,3) 

where : 


a,  +  a      -    1  (C.4) 


and: 


0  _<  a,   <   1 

0  <  a2      <       1  (C.5) 

Equation  C.3  reduces  to  the  small  deflection  profile,  C.l,  for 
a.  =  1,  a?  =  0;  and  to  the  membrane  profile,  C.2,  for  a,  =  0, 
a.  =  1.   Using  this  profile  and  proceeding  as  in  Chapter  II: 
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4      3r2         2 
K   =  -w"  =  w  [a      ±r    (l-^)  +  a_  ^-]  (C.6) 

r  o   1  a      a        2  a 

KQ  =  -—   =  w  [a,  ^(1-^)  +  a_  ^-]  (C.7) 

0     r      o   1  a    a      2a 


and 


W  2 

Kr  +  ±<q   =  5?  [2a1(3-7|r)  +  3a2]  (C.8) 

W  2 

ka  +  *<„  =  -§  [2a,(3-5§x)  +  3a9]  .  (C.9) 

y      r    a     x     a       z 

In  these  equations  and  those  to  follow  we  may,  of  course, 
eliminate  either  a.  or  a.  by  using  C.4.   This  was  not  done 
here,  however,  in  order  to  keep  the  terms  from  the  C.l  profile 
and  those  from  the  C.2  profile  directly  visible  and  separated 
in  each  equation. 

For  the  radial  deflection  rate,  u,  we  use  the  form: 

u  =  r(a-r)  (c,+c2r)  (CIO) 

The  extensional  deflection  rates  are: 

e   =  u'  +  w'w'        eQ  =  -I  (C.ll) 

r  0    r 

Using  as  before: 

WW  WW* 

r        ,   o  o         ,   o  o 
X  =  a'  cl  °  ci  T*-  '  C2  "  C2  "T5-  < 

w  w 

e      =  _P_^>  f(i-2x)c'  +  (2x-3x2)ci  +  16x2a2  (1-x2 )  2 
r     a  1  2        1 

+  16x2a1a2(l-x2)  +  4x2a2]  (C.12) 

w  w 
ee   =  -|^  [■•(i-x)cj-  +  (x-x2)c2]  (C.13) 
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These  reduce  to  the  appropriate  expressions  for  the  small 

deflection  profile  when  a.  =  1,  a_  =  0,  and  to  the  appropriate 

membrane  profile  expressions  for  a,  =  0,  and  a_  =  1. 

w  w 
*r  +  *«8  =  "f^2  [(l  "  *x)cl  +  <*x-***c2  +  16x2aJ(l-x2)2 


+  16x2a1a2(l-x2)  +  4x2a2] 


(C.14) 


w  w 
+  fer    -  ~A    [(f  -  2x)c^  +  (2x  -  |x2)c^ 


+  8x2a2(l-x2)2  +  8x2a1a2(l-x2)  +  2x2a2](C15) 


The  general  expression  for  moment  per  unit  length  derived 

earlier  is :  ; 


„     ,4  ,1/n   n    1  rh.  ,bhil+l/n        n 
M  =   {3K}     n+T  b  [2"la+2-|       "  (2n+l)b 


bh,2+l/n 
a+2"'s 


hj  _bhil+l/n      n    ,    bh.2+l/n, 
+  2|a  2~'       +  (2n+l)b'  a  2~ Is     J 


(C.16) 


and  for  extensional  forces  is: 

n    1  ri  .bhil+l/n    i   bh  1 1+1/n,    ,_  ,  _,v 
5+1   b  (la+2-|       -  la"2-|      ]  <C-17) 


N=  4,1/" 


Where  for  M   and  N 
r      r 


a  =  er  +  7efl 

b  -  K   +  iftj 

r    ^9 


and  for  M   and  N 


a  =  efl  +  7e 

b  =  k Q  +  £k 
0      r 


Defining: 


A   =  2a, (3  -  7X2)  +  3an 
c      1  2 
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w  w 


Bc   =    (I   -    §x)c{   g2     +    (|x  -    ^x2)c'      ^ 


w  w 

+    16x2a2(l-x2)2    Wo      +16x2a.,a.,(l-x2)    ^  +    4x2a2    r^ 
1                    r—                    12                  h  2   h 

n 


D      =    2an  (3   -    5x2)    +    3an 

C  1  2 


WW  w 

Ec    =     <*    -    2x)c'    ^  +    (2x    -    fx2)c'    g2      +    8x2a2(l-x2)2    ^ 


V7  W 


+    8x2a,a0(l-x2)    r-£  +    2x2a2    r^    , 
12  h  2   h 


we   may  write: 


„   J      =    (4    ,1/n      n      h2+l/n      "o    1+1/n  1       ,  j       ( l+l/„ 

r  r     3K      n+1  a  2A   '  c      c' 

c 


n       i_     i-  i2+l/n     1  |  „     1.  i 1+1/n 

(2n+l)  A*   lBc  +  ^Acls       +  2A-I  Bc  "  *Ac  I 
c  c 


2+1/n,  .  ry.   ,,  ,v2 


+(2n+l)  A^  I  Bc  "  *Acls     !  '  f4al(l-3x2)+2a2](C.18) 


N  4  =  (l^)1/"  -a,.  h2+1/"  (!°)1+1/n  i  {|b  ♦  iA  |1+1/n 

r  r     3K      n+1  a2        A    '  c      c1 

c 

w  w 

-|B   -  iAc|1+1/n]  '  t(l-2x)c«  ^  +  (2x  -  3x2)c'  ^ 


w  w 

+  16x2«|(l-x2)2  ^  +  16a1a2(l-x2)x2  ~ 


w 
+  4x2a|  — ]  (C.19) 
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•   _  ,4_.l/n  _n_     ,2+1/n  ,^o.  1+1/n  ,  1  ■  _,     in  .  1+1/n 


n     i  „      A  in  i2+l/n  .  _1  i  «     i  r,  1 1+1/n 


E 


(2n+l)D^  '   c    '  c's        2D  ■   c 
c  c 


I  t,     iD  1 2+l/n]  .  [4a  (1_x2}  +  2an](C20) 


+  (2n+l)Dz  I  Ec    * 
c 


Ve  =  <fe'1/n  £l    h2+Vn  <^»1+1/n  o^  Mbc  +  *Dct™ 

C 

14-1/  W  W 

-|Ec  -  iDc|1  +  1/n]  -  [(l-x)c^  ^   +  (x-x2)c!>  ^].(C21) 


The  total  energy  dissipation  rate  is  obtained  by  integrating 
these  expressions  over  the  deformed  surface  area  of  the  plate. 
This  gives: 

V  =  27ra2  (|^)1/n  nTI  ^^    (^?)  1+1/n  Ic  (C.22) 


where  I   is  the  integral  of  the  bracketed  expressions  above 

w 
over  x  for  a  particular  case  (r— -*  n) ,  and  where  the  integral 

has  been  minimized  with  respect  to  c^  and  c'  for  the  value  of 

a.  and  a?  under  consideration. 

Takinq  a  small  variation  of  w   about  its  equlibrium 

o 

w 

position  for  some  particular  value  of  j — : 

r{\7  •  o     . 

££-      6w        =      2tt    /q    6w   dr  (C.23) 

d\\r  o  6  r 

o 


or 


•  1/n  r~   2,4  ,1/n  ,2+1/n  ,1  , 1+1/n  _  ,  r • 

[2na  (3K}     h       (a^>       V  6Wo  = 


*2  +  cu(l-?4-)]rdr.  (C.24) 


2nq6w^   /a  [a,  (1-rr) 2  +  a0  (I-tt)  ]  rdr 
o   o     -L    a       z    a 
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Solving  for  w  : 

2+2/n 

*o   =  [  ,4  ,l%\2;i/n   12 1 ]U    >  <C'25> 

(3K}     h       3=^7  Xc 

or 

f,a.2+l/n  ,3Ka*l/n    ,    al  ,,n  ,„  „rs 

o   =  f  (h)       (~1~)     q  (12~I — )J   '  (C.26) 

c 

and  in  non-dimensional  form: 

w  0      -,  ,  i    3-a, 

(   °  )  =  1  (a}2n+l  ( l_}n                       ,     . 

V   n;  4  V      ^12  I   '                        Kl*.*f) 

Kaq  c 
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APPENDIX  D 
Mathematical  Details  for  Analysis  of  a  Simply  Supported 

Circular  Plate 

In  calculating  the  elastic  deflection  of  a  simply 
supported  circular  plate  subjected  to  a  uniform  lateral 
pressure  Timoshenko  and  Woinowsky-Krieger  [26]  find  that  the 
deformation  profile  assumes  the  following  form: 

w  =  wn    (1-  K)   (1  -  III      §4)  (D.D 

o      a        5  +  v   a 

where  w   is  the  center  deflection.   If  we  make  v  =  1/2  and 
o 

take  derivatives  with  respect  to  time,  we  get  the  deformation 
rate  profile: 

*  «  *„  <1  - '  fr>  (1  -  TT  §*>  •  (D-  2) 

o      a       11   a 

One  may  now  compute : 

.  „     o  r28    36   r2  t  rr\    ->\ 

Kr  =  "W   =  I*  [TT  "  IT   a"^1  (D*3) 

• 

w'    Wo  r28    12   r2,  ,n  A, 

Ke  =  "P-  =  iT  fn  "  II  I*"!'  (D*4) 


and 


i .      o  r 42    42  r2 ,  ,n  n 

Kr  +  N  =  ^    [IT  -  IT  a^]  (D*5) 

J  i .     Wo  r42    30  r2,  ,n  Cx 

K9  +  7Kr  =  P"  [TT  -  IT  I**-  (D'6) 


Using,  as  in  Chapter  II,  the  expression: 

u  =  r(a-r)  (cJL    +   c2r)  (D.7) 


for  radial  deflection  rate  where, 
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WW  WW 

cl  -  ci  "P2  C2  -  C2  -f^  (D-8» 

and  letting  x  =  r/a: 

w  w 

e   =  u'  +  w'w'  =  -^  [(l-2x)c*  +  (2x-3x2)c: 
r  a  1  2 

+  x2(ff  "  j~   x2)2]  (D.9) 

W  W 

e0  =  £  =  -^r-      [(l-x)c^  +  (x-x2)c']  (D.10) 


and 

6r  +  *6e  =  "IT^  [(*   *x)ci  +  (*x   ?x2>c2 


w  w 
o  o 


2 /28     12   2 \ 2 


+x z(ff  "  ffxz)z]  (D.ll) 


W  W 

ee  +  *er  =  -|^  [(1  -  2x)c^  +  (2x  -  fx  )  c£ 


+  ~  (ff  -  YT  x2)2]  (D*12) 


These  expressions  may  now  be  substituted  into  the  equations  for 
moment  and  extensional  force  resultants,  (2.14)  through  (2.17). 

Defining : 

A   =  (il  _  42   2) 
s  ~  lll    11  x  ' 

WW  w 

Bs  -    (*       7X)iTci  +   (7X      7X  )iTc2  +  x   (n       nx  ]  h 

42         30    __2 


Ds  E    (IT-  ITX  > 


W  W  9         ««  t  ^  w 


Bb    3    .(«   -    2.x,c1    ^  *    (2,    -    IxMc'    ^  +  5-    (fi  -    «««,.^ 
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we  may  write: 

w 

..  •     /4  x  1/n   n   ,2+1/n  .  o* 1+1/n  r  1  i  _     i  _  i 1+1/n 

M  <       =  (^77)      — TT    h        (~r)  [o5T~  B   +  ?A       / 

r  r     3K      n+1          a  2A  '  s      s1 


n     i„     1.  i2+l/n     _1  I  R     1,  I 1+1/n 

-7-5; — .  t  \  ■»  2     B  +    7  A  +    tt^ —  B       -     7  A 
(2n+l)A    '    s              s's                      2A    '       s  s1 

s  s 


,  n iD  1,     i2+l/n,  r28        36    „2i  /~    n\ 

+   -7-^ — ,  t  .  .;  2    B      -    7A  J     *     Itt"  ~    tt"  x    J  (D.13) 

(2n+l)A'£  '    s         z    s's  11         11 


„   e      -    (4=)    -5r  h2+1/"    (V+1/n  -i-    [|B  +   iA    i1+1/n 

r   r           3K      n+1                        a^                     A         '    s  *    s1 

s 

-|B3    -    iAs|1+1/n]     •[(!'-    2x)c'    ^  +  (2x   -    3x2)c«    ^ 


o 


>  ,  28    12  2  v  2   o 


+  xMfy  -  ff*2)'  ^]  (D.14) 


e 


*   -  (*-)       n   h2+l/n  ,  o. 1+1/n    1  I      1   1 1+1/n 

K6  "  (3K)  nTI  h       (a^r)       [2D~  |Es  +  7°s  ' 

s 

^^2    |e     +  Id   |2+1/n  +  ^-  |e     -  io   |1+1/n 

!n+l)D^     '    s         *    s's  2D  '    s  s1 


(2n+L. 

s  s 


+  TzhW  lEs  -  ±°sls+VnJ  •  [ff-Tr-2]         <D-15' 


s 

,  . . ' .  w  w 

-|ES    -    iDg|1  +  1/n]     •     [(l-x)c'       £      +    (x-x2)^^]. 

(D.16) 
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The  internal  energy  dissipation  rate  is: 


V  =  27ra2  /l(MvKr  +  MflKft  +  N  £   +  N_efl)  xdx  (D.17) 

o    r  r    do    rr    u  u 


giving  with  (D.13)  through  (D.16): 

v        -5^^  2  /4  ,1/n  n   .  2+1/n  ,  o.  1+1/n  T  ,n  ,  Q> 

V  =  2lTa   (3K)     nTT  h       (iPr)       Xs  (D*18) 


where  I   is  the  integral  of  the  bracketed  geometric  terms  in 
the  above  expression.   From  the  principle  of  virtual  velocities 


—  5w   =  2tt   /aq  Sw  rdr  (D.19) 

j  •     o        o 

dw 
o 

which  gives: 

.  1/n  ,1  . 1+1/n   „   2/4  ,1/n   ,  2+1/n  _   x  . 
w      (—7-)    /    2?Ta  (^7)     h    /   I   6w   = 
o     a^  3K  s    o 


27rq  6w   /a  (1-^2-)  (1~tt  ^t)  rdr  =  2irq6w   ~  a2  .   (D.20) 
o  0     a     11  a  o  2.  a 

So: 

*o  ■  '4  ^lill/n   22  T  »"  (D-21> 

<3k'     h      5-  *s 


,-,      _    r /a.  2+1/n    .3Ka.  1/n  1         ,n  ,      00v 

or        WQ=[(R)  (~r)  q    2_2  1  (D.22) 

5        s 


°r  ~^=  <h->2n+1  I1277— '"•  (D-23) 

Kaq  -= —  I 

5        s 
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APPENDIX  F, 
Illustration  of  Technique  for  Finding  the  Minimum  Energy  Integral 

As  explained  in  Chapter  II,  the  quantities  c'  and  cl    are 

determined  by  requiring  that  the  internal  energy  dissipation 

rate  (2.37)  be  minimized.   Since  the  integration  required  to 

evaluate  (2.37)  must  be  done  numerically,  this  minimization 

must  be  accomplished  through  a  systematic  search.   The 

w 
minimization  must  be  done  for  each  case  (zr~t    n)  of  interest; 

but,  by  considering  a  systematically  varied  series  of  cases, 

a  set  of  curves  such  as  those  in  Figures  2.2  and  2.4  can  be 

obtained. 

The  procedure  used  for  minimization  can  be  illustrated 

w 
by  considering  a  typical  case,  n  =*  5.0  and  r—  =  1.0  for 

clamped  edges.   The  computer  program  listed  in  Appendix  H, 

when  provided  an  initial  value  for  c^  and  a  value  for  c' , 

will  compute  the  integral,  I  ,  for  the  given  c,1  and  c'  values 

at  regular  intervals  of  0.3.   In  the  accompanying  plots,  the 

quantity  $  rather  than  the  integral  I   is  plotted  so  that  $ 

may  be  transferred  directly  to  Figure  2.2.   As  $  =  f(j~ )>  the 

c 
maximum  value  of  ^  corresponds  to  the  minimum  value  of  I  . 

As  seen  in  the  plots,  c'  is  varied  until  a  peak  is  found  for 

a  particular  value  of  c'.   Peaks  are  then  found  for  other 

values  of  ci  until  a  maximum  is  found  with  respect  to  c'.   The 

resulting  envelope  of  curves  gives  a  maximum  value  of  $  which 

gives  a  point  on  the  plot  of  Figure  2.2.   The  procedure  is 

w 
then  repeated  for  the  other  values  of  r—  and  n. 
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There  is  no  mathematical  assurance  that  the  value  of  $ 
found  as  above  represents  a  unique  and  absolute  maximum  value. 
The  initial  area  of  search  was  that  suggested  for  the  analytical 
solution  of  the  n  =  1.0  case;  and  a  very  broad  parameter  search 
with  respect  to  c,'  and  ci  gave  no  evidence  of  any  higher 
maxima.   The  agreement  between  the  experimental  results  and 
the  theoretical  results  suggests  that  the  values  of  $  obtained 
are,  in  fact,  the  proper  ones. 
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APPENDIX  F 
Tabulation  of  Data  from  Plate  Experiments 

The  following  tables  list  the  data  taken  during  the  plate 
experiments  discussed  in  Chapter  I. A.   All  tests  were  taken  at 
room  temperature  which  generally  was  75°F  ±3°F  with  some  points 
outside  these  limits.   It  is,  of  course,  a  well  known  fact  that 
creep  behavior  in  metals  is  affected  by  temperature.   Within 
the  relatively  narrow  temperature  band  here  considered,  however, 
such  effects  may  be  expected  to  be  rather  small  and  definitely 
secondary  to  the  influence  of  pressure  and  geometry  effects 
which  were  to  be  investigated  in  the  plate  tests.   Examination 
of  the  data  reveals  no  appreciable  systematic  influence  of 
temperature  in  these  tests.   All  times  are  in  hours  and  minutes, 
all  deflections  are  in  inches,  all  pressures  in  pounds  per 
square  inch,  and  all  temperatures  in  °F. 
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TEST 

1 

8  LB. 

PLATE 

3, 

.  0  psi 

Time 

Temp 

Hrs-Mins 

r=0 

r=4" 

r=l" 

r=l£M 

r=2" 

r=2*" 

°F 

00 

.0052 

72 

01 

53 

.0049 

.0033 

.0021 

.0005 

.0001 

03 

54 

49 

33^ 

21 

05 

0l£ 

15 

57 

53 

36 

22 

06 

oij 

32 

58* 

55 

37 

24 

06* 

oii 

72 

50 

60 

57 

38* 

24 

07 

02 

1-00 

61 

57 

39 

24 

07 

02 

72 

1-30 

63 

59 

40 

25 

07^ 

02 

2-00 

64 

61 

41* 

26 

09 

02 

72 

2-42 

66 

63 

43 

27 

09 

02 

3-00 

67 

64 

43 

27^ 

09 

02 

72 

4-45 

71 

67 

61 

29 

10 

02 

5-20 

72$ 

68 

47 

29| 

10 

02* 

72 

6-00 

73 

68$ 

48 

30 

10 

02} 

6-30 

74 

70 

48^ 

31 

10 

02j 

72 

25-00 

92 

88 

62 

39 

14 

07 

74 

47-00 

.0110 

.0105 

75 

47 

17 

08 

76 

52-00 

14 

08 

78 

48 

18 

08* 

76 

73-00 

23 

20 

86 

53 

20 

10 

76 

100-30 

30 

26^ 

90 

55* 

21 

10 

76 

125-40 

34 

30 

93 

57 

21* 

11 

76 

TEST  2 
8  LB.  PLATE  3.0  psi 

Time       r=0    r=4"    r=l"    r=l£"    r=2"    r=2*"    Temp 


00  .0052 

01  53   .0045   .0036   .0020    .0007   .0002      77 


03 

54 

45 

37 

21 

07 

02 

15 

56 

51 

38* 

23 

08* 

02 

32 

58 

53 

40 

23* 

09 

02 

76 

50 

59 

53 

41 

24 

09 

02 

1-00 

59^ 

53i 

41 

24 

09 

02 

76 

3-50 

65 

62 

45 

28 

11 

02 

76 

4-30 

67 

64 

46 

28 

Hi 

02| 

76 

5-00 

68 

64$ 

47 

29 

12 

03 

76 

5-32 

69 

65 

48 

29 

12 

03 

76 
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TEST    3 


8  LB 

.  PLATE 

4.1  psi 

Time 

Temp 

Hrs-Mins 

r=0 

r=p 

r=l" 

r=l£ 

r=2" 

r=2fc" 

op 

00 

.0352 

01 

.0388 

.0353 

.0280 

.0181 

.0060 

.0023 

69 

03 

.0405 

70 

93 

89 

63 

25 

07 

26 

85 

.0304 

97 

65 

26 

11 

33 

94 

11 

.0201 

67 

26 

18 

44 

.0404 

19 

05 

69 

28 

30 

56 

15 

28 

10 

71 

29 

68* 

47 

68 

26 

36 

16 

73 

29 

1-00 

75 

33 

41 

19 

74 

30 

68£ 

1-26 

84 

41 

47 

23 

75 

31 

1-49 

93 

48 

54 

27 

77 

31 

69 

2-30 

.0510 

64 

66 

34 

80 

36 

3-00 

17 

70 

72 

37 

81 

36 

69 

4-10 

32 

84 

81 

44 

83 

37 

5-00 

40 

91 

88 

48 

85 

37 

69 

5-30 

46 

98 

92 

51 

86 

.   38 

5-50 

49 

.0500 

94 

52 

87 

38 

69£ 

TEST    4 


8  LB. 

PLATE 

4.1  psi 

me 

r=0 

r=r  ; 

r=l" 

r=liT 

r=2" 

r=2£" 

Temp 
°F 

00 

.0310 

70 

01 

.0355 

.0311 

.0251 

.0145 

.0055 

.0017 

03 

63 

25 

62 

51 

58 

18 

06 

78 

38 

77 

57 

60 

19 

11 

91 

50 

83 

61 

62 

21 

15 

98 

56 

88 

64 

63 

21 

70 

22 

.0407 

64 

94 

67 

64 

21 

30 

18 

73 

.0300 

72 

65 

22 

45 

30 

85 

09 

77 

68 

22 

1-00 

37 

91 

15 

80 

69 

23 

70 

1-41 

56 

.0410 

28 

88 

72 

24 

2-00 

62 

15 

32 

91 

74 

25 

70 

2-30 

71 

24 

39 

95 

76 

25 

3-00 

79 

31 

45 

97 

77 

25 

70£ 

3-33 

87 

39 

50 

.0201 

78 

25 

4-06 

94 

45 

56 

04 

80 

26 

71 

4-30 

.0500 

50 

60 

07 

80 

27 

5-00 

05 

55 

64 

09 

82 

27 

71 
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TEST 

5 

8  LB. 

PLATE 

6.1  p: 

si 

Time 

Hrs-Mins 

r=0 

r=4" 

r=l" 

r=lj-r 

r=2" 

r=2$- 

00 

.0824 

01 

.0869 

.0805 

.0614 

.0376 

.0150 

.0055 

03 

96 

28 

32 

87 

55 

57 

06 

.0916 

47 

46 

95 

59 

59 

15 

44 

75 

66 

.0408 

64 

62 

25 

61 

94 

81 

17 

68 

64 

31 

68 

.0902 

86 

21 

70 

64 

45 

87 

15 

98 

27 

74 

66 

1-00 

99 

28 

.0708 

34 

77 

66 

1-30 

.1020 

48 

23 

44 

80 

68 

2-00 

35 

64 

34 

52 

84 

70 

3-00 

59 

89 

53 

64 

90 

72 

3-36 

72 

.1000 

62 

69 

93 

73 

4-45 

93 

20 

76 

79 

98 

75 

5-15 

.1107 

27 

86 

83 

.0200 

82 

6-00 

19 

40 

93 

88 

03 

82 

Temp 


79 


79 


80 


80 


81 


81 


TEST 

6 

8  LB. 

PLATE 

6.1  psi 

me 

r=0 

r=£" 

r=l" 

r=l£" 

r=2" 

r=2j-" 

Temp 
°F 

00 

.0840 

01 

.0883 

.0842 

.0631 

.0379 

.0156 

.0073 

76£- 

03 

.0908 

65 

49 

90 

61 

75 

07 

33 

86 

66 

.0400 

65 

78 

15 

55 

.0910 

83 

11 

70 

80 

34 

85 

36 

.0707 

24 

76 

85 

76^- 

48 

98 

49 

17 

30 

79 

86 

1-00 

.1006 

57 

22 

34 

80 

87 

76^- 

1-25 

21 

70 

35 

41 

83 

89 

1-50 

33 

81 

44 

47 

87 

92 

77 

2-00 

38 

87 

48 

49 

88 

93 

3-15 

61 

.1009 

68 

59 

92 

95 

78 

4-20 

78 

23 

81 

66 

96 

97 

5-50 

98 

43 

98 

76 

.0201 

.0100 

78 

106 


TEST 

7 

8  LB. 

PLATE 

8.01  psi 

Time 

Temp 

Hrs-Mins 

r=0 

, l  it 

r=l" 

r=li" 

r=2" 

r=2£" 

°F 

00 

.1270 

73 

01 

.1313 

.1251 

.0943 

.0613 

.0244 

.0120 

03 

39 

76 

65 

28 

50 

24 

06 

65 

98 

82 

39 

55 

29 

12 

88 

.1321 

.1002 

52 

60 

32 

73 

34 

.1429 

57 

36 

73 

70 

38 

1-00 

54 

84 

57 

88 

76 

42 

1-27 

75 

.1403 

74 

98 

81 

46 

73 

2-00 

91 

22 

88 

.0708 

86 

48 

3-23 

.1524 

49 

.1114 

26 

93 

54 

73 

4-00 

34 

58 

22 

32 

96 

55 

5-15 

52 

76 

36 

42 

.0301 

58 

74 

TEST 

8 

8  LB. 

PLATE 

8. 

.01 

psi 

me 

r=0 

r=i" 

r=l" 

r=li" 

r= 

=  2" 

r=2V 

Temp 

op 

00 

.1252 

76 

01 

.1301 

.1247 

.0921 

.0586 

.02: 

.0095 

03 

28 

74 

42 

99 

43 

98 

05 

43 

88 

55 

.0607 

47 

99 

18 

89 

.1333 

91 

33 

57 

.0104 

30 

.1413 

54 

.1012 

45 

62 

10 

76 

45 

31 

70 

28 

55 

66 

12 

1-00 

45 

83 

39 

63 

69 

14 

76 

1-15 

56 

93 

48 

69 

72 

15 

1-30 

65 

.1403 

56 

75 

75 

16 

2-20 

89 

23 

76 

87 

80 

19 

76 

3-06 

.1504 

40 

90 

97 

85 

21 

3-35 

17 

51 

98 

.0703 

88 

24 

75 

4-00 

24 

57 

.1104 

07 

89 

24 

76 

4-45 

36 

68 

14 

15 

93 

27 

76 

5-45 

50 

83 

26 

23 

96 

30 

6-00 

53 

86 

29 

25 

97 

30 

76 
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TEST    9 
6    LB.    PLATE  2.39    psi 


Time 

Temp 

Hrs-Mins 

r=0 

r=!" 

r=l" 

r=l!" 

r=2" 

r=24" 

°F 

00 

.0099 

74 

01 

.0100 

.0091 

.0067 

.0025 

.0012 

.0001 

09 

12| 

.0102^ 

74^ 

31 

134 

oil 

15 

15* 

05 

77 

32 

14 

oil 

21 

18 

07! 

78$ 

33 

144 

02 

30 

20 

11 

80 

344 

15 

02 

45 

23 

13* 

82| 

36 

16 

03 

74 

1-00 

25| 

14 

84 

374 

16| 

03 

1-15 

274 

17 

86 

38 

164 

03 

1-30 

294 

19$ 

874 

384 

17 

03! 

2-00 

32^ 

23-1 

90 

404 

18 

05 

74 

3-17 

37 

30 

93i 

43 

19 

05 

4-00 

384 

3l| 

95 

444 

19 

054 

74 

4-30 

39! 

32 

96^ 

444 

19 

054 

4-45 

394 

32 

96! 

444 

19 

05! 

5-00 

41* 

334 

97 

454 

194 

05l 

74 

5-15 

4l4 

334 

97! 

45l 

194 

05! 

5-30 

41-S 

33| 

974 

45! 

19! 

05! 

6-15 

4l| 

364 

99 

47 

20 

05! 

74 

TEST 

10 

6  LB 

.  PLATE 

2. 

.39  psi 

me 

r=0 

r=l" 

r=l" 

r=l!" 

r=2" 

r=24" 

Temp 
°F 

00 

.0097 

74 

01 

.0100 

.0089   , 

.0068 

.0026 

.0010 

.0000 

09 

09 

.0100 

75 

31 

12 

00 

74 

15 

12 

03 

77! 

32 

12 

oo! 

21 

14 

05 

79 

33 

13 

01 

30 

16! 

07 

80l 
82* 

34 

13 

01 

46 

19! 

11 

36 

13! 

01 

74 

1-00 

21! 

13 

84 

37 

14 

01 

1-15 

23 

13! 

85 

38 

14! 

01 

1-30 

244 

144 

86 

38 

144 

01 

2-00 

271 

18 

88f 

924 

394 

15 

02 

74 

3-15 

33 

23 

42 

16 

024 

4-00 

34! 

254 

94 

43! 

17 

02! 

4-30 

35 

27 

944 

44 

17 

024 

5-00 

354 

28 

95! 

444 

17 

02* 

5-15 

354 

28 

954 

44^ 

17 

02* 

5-30 

351 

284 

96 

45 

17 

03 

6-15 

354 

294 

97 

454 

17 

03 

108 


TEST  11 
6  LB.  PLATE  4.1  psi 


Time 

Temp 

Hrs-Mins 

r=0 

r=*" 

r=li" 

r=l|M 

r=2" 

°F 

00 

.0510 

76 

01 

57 

.0498 

.0350 

.0265 

.0081 

06 

79 

27 

71 

80 

86 

12 

94 

42 

81 

87 

89 

15 

.0600 

47 

85 

90 

90 

76 

21 

09 

54 

91 

94 

92 

30 

19 

63 

98 

99 

94 

76 

45 

31 

74 

07 

05 

96 

1-00 

40 

82 

13 

09 

97* 

76 

1-15 

47* 

88 

18 

13 

00 

1-30 

52 

93 

22 

15 

00* 

2-00 

62 

03 

29 

20 

02* 

76 

2-30 

71 

11 

35 

24 

04 

3-00 

79 

20 

41 

28 

06 

76 

4-02 

90* 

32 

49 

34 

08* 

4-30 

97£ 

36 

53 

37 

09 

76 

4-45 

99 

39 

55 

38 

10 

5-00 

.0701 

41 

57 

40 

12 

5-15       03*     43      59      41*      12     76 


TEST  12 
6  LB.  PLATE  4.1  psi 


Time 

Temp 

Hrs-Mins 

r=0 

r=*" 

r=l" 

r=l*" 

r=2" 

°F 

00 

.0530 

76 

01 

.0574 

.0519 

.0399 

.0247 

.0092 

06 

12 

52 

24 

63 

98 

12 

26 

67 

34 

71 

01 

76 

15 

33 

71 

37 

73 

02 

24 

43 

81 

44 

79 

05 

30 

52 

89 

50 

82 

06 

53 

72 

.0607 

65 

92 

10 

76 

2-22 

.0708 

41 

94 

.0310 

17 

76 

2-45 

15 

47 

99 

13 

18 

3-00 

21 

51 

.0502 

15 

19 

76 

3-31 

29 

59 

09 

20 

21 

4-30 

40* 

70 

18 

26 

23 

5-00 

44 

75* 

23 

28 

24 

76 

5-30 

51 

79* 

27 

31 

25 

6-00 

54 

83 

30 

34 

26 

76 

6-12 

57 

85 

31 

34 

26* 
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TEST  13 
6  LB.  PLATE  6.1  psi 


Time 

Temp 

Hrs-Mins 

r=0 

r=£" 

r=l" 

r=l£" 

r=2" 

r=2*" 

op 

00 

.0880 

76 

01 

.0927 

.0852 

.0650 

.0441 

.0175 

.0071 

03 

51 

74 

68 

53 

80 

73 

06 

68 

90 

80 

62 

84 

75 

11 

85 

.0905 

94 

71 

88 

78 

15 

94 

15 

.0701 

76 

90 

79 

30 

.1017 

38 

20 

90 

96 

81 

76 

40 

29 

49 

28 

96 

98 

83 

45 

33 

54 

33 

99 

.0200 

83 

1-00 

45 

64 

42 

.0505 

03 

85 

76 

1-17 

57 

78 

53 

11 

05 

87 

1-31 

66 

84 

59 

15 

06 

89 

77 

2-00 

83 

.0100 

71 

24 

11 

91 

2-30 

94 

11 

80 

31 

13 

92 

77 

3-00 

.1104 

22 

89 

37 

15 

93 

ll\ 

4-15 

26 

43 

.0807 

50 

21 

95 

4-45 

33 

51 

14 

54 

23 

96 

ll\ 

5-30 

45 

60 

23 

60 

26 

.0100 

TEST  14 
6  LB.  PLATE  6.1  psi 

Time       r=0    r=\"        r=l"    r=l7H    r=2"    r=2TH   Temp 

op 


80 

80 
80 

80 

80 
80 


00 

.0878 

02 

.0936 

.0865 

.0643 

.0438 

.0152 

.0078 

06 

62 

88 

65 

51 

57 

82 

15 

87 

.0915 

85 

65 

62 

86 

31 

.1016 

42 

.0706 

79 

68 

97 

45 

32 

56 

18 

88 

71 

98 

1-00 

44 

68 

28 

94 

74 

.0100 

1-23 

60 

83 

40 

.0503 

78 

03 

1-30 

66 

88 

45 

06 

79 

04 

2-00 

79 

98 

55 

13 

81 

05 

2-30 

91 

.1013 

66 

21 

84 

07 

2-45 

97 

17 

70 

24 

84 

07 

3-16 

.1106 

26 

77 

30 

86 

08 

4-00 

18 

37 

87 

36 

89 

09 

5-00 

30 

48 

98 

44 

92 

11 

5-30 

36 

56 

.0802 

48 

94 

12 

110 


TEST 

15 

6  LB. 

PLATE 

8.0  psi 

Time 

Temp 

Hrs-Mins 

r=0 

r=4" 

r=l" 

r=li" 

r=2" 

r=2£" 

°F 

00 

.1308 

75 

01 

.1354 

.1265 

.1009 

.0679 

.0255 

.0114 

03 

78 

88 

28 

93 

61 

16 

06 

98 

.1308 

44 

.0705 

66 

19 

15 

.1427 

36 

72 

25 

74 

23 

75 

24 

57 

54 

84 

36 

80 

27 

30 

63 

61 

92 

42 

83 

29 

75 

45 

73 

78 

.1108 

54 

88 

31 

1-00 

87 

91 

19 

63 

92 

33 

75 

1-15 

97 

.1403 

29 

70 

96 

35 

1-36 

.1510 

14 

40 

79 

.0300 

36 

75* 

2-00 

24 

27 

51 

88 

05 

38 

2-45 

44 

51 

69 

.0803 

11 

44 

76 

3-00 

50 

59 

74 

07 

12 

44 

4-20 

74 

81 

97 

25 

20 

47 

76 

5-04 

87 

95 

.1206 

33 

24 

48 

5-30 

93 

.1500 

12 

38 

26 

49 

76? 

6-00 

.1600 

06 

18 

43 

28 

50 

TEST    16 


6  LB. 

PLATE 

8.0  psi 

me 

r=0 

r=£" 

r=l" 

r=l£" 

r=2" 

r=2i" 

Temp 
°F 

00 

.1295 

76 

01 

.1342 

.1248 

.1009 

.0636 

.0265 

.0087 

03 

68 

78 

30 

53 

72 

89 

06 

88 

97 

47 

65 

77 

92 

15 

.1417 

.1322 

73 

83 

86 

96 

76 

24 

34 

44 

88 

97 

91 

98 

45 

65 

68 

.1113 

.0716 

.0300 

.0103 

1-00 

78 

85 

25 

27 

05 

04 

76 

1-16 

90 

98 

37 

35 

09 

06 

1-30 

99 

.1406 

44 

42 

12 

07 

2-00 

.1517 

19 

60 

53 

18 

09 

76 

2-30 

31 

36 

71 

63 

23 

10 

3-04 

44 

47 

83 

73 

27 

12 

3-30 

53 

59 

92 

80 

30 

13 

76 

111 


TEST 

17 

4  LB. 

PLATE 

0.82  psi 

Time 

Temp 

Hrs-Mins 

r=0 

r=i" 

r=l" 

r=li" 

r=2" 

r=2|" 

°F 

00 

.0110 

75 

01 

.0112 

.0106 

.0080 

.0050 

.0018 

.0003 

03 

15 

08 

82 

50 

18 

03 

06 

18 

13 

84 

52 

20 

04 

15 

23 

16 

88 

53 

20 

04 

741 

30 

27 

21 

91 

56 

21 

04 

1-00 

33 

26 

94 

58 

22 

04 

75 

1-34 

39 

32 

98 

59 

23 

05 

3-00 

45 

38 

.0103 

63 

25 

05 

74l 

4-15 

53 

46 

08 

66 

26 

06 

5-04 

56 

48 

10 

67 

27 

06 

74! 

6-00 

57 

50 

11 

68 

27 

06 

TEST 

18 

4  LB, 

.  PLATE 

0.82  psi 

me 

r=0 

r=£" 

r=l" 

r=li" 

r=2" 

r=2£" 

Temp 
°F 

00 

.0127 

73 

01 

.0131 

.0122 

.0080 

.0055 

.0020 

.0018 

03 

33 

25 

81 

56 

21 

19 

06 

35 

27 

83 

57 

23 

20 

15 

39 

31 

86 

59 

22 

20 

73 

30 

42 

33 

87 

60 

23 

21 

1-00 

45 

37 

89 

61 

24 

22 

73 

1-30 

47 

39 

89 

62 

24 

22 

2-00 

49 

42 

91 

63 

25 

22 

2-18 

49 

42 

91 

63 

25 

22 

73 

112 


TEST 

19 

4  LB 

.  PLATE 

1.840  p; 

si 

Time 

Temp 

Hrs-Mins 

r=0 

r=£" 

r=l" 

r=l£" 

r=2H 

r=2fc" 

°F 

00 

.0845 

76^- 

01 

.0871 

.0820 

.0653 

.0452 

.0170 

.0045 

03 

83 

32 

64 

55 

73 

46 

06 

92 

41 

72 

61 

75 

47 

15 

.0906 

55 

83 

68 

79 

48 

21 

11 

60 

88 

71 

81 

48 

31 

18 

68 

93 

76 

83 

49 

76^ 

45 

24 

74 

98 

80 

84 

49 

1-00 

29 

80 

.0702 

83 

86 

50 

76£ 

1-31 

34 

86 

08 

88 

87 

50 

2-00 

42 

92 

13 

92 

89 

51 

76* 

2-30 

47 

98 

18 

95 

91 

51 

3-00 

51 

.0901 

22 

98 

92 

52 

77 

TEST    20 


4  LB 

.  PLATE 

1.840  p; 

si 

me 

r=0 

r=£" 

r=l" 

r=li" 

r=2" 

r=2i" 

Temp 
°F 

00 

.0820 

79 

01 

.0848 

.0783 

.0638 

.0438 

.0162 

.0053 

03 

61 

98 

49 

46 

65 

55 

06 

70 

.0804 

56 

50 

68 

56 

12 

79 

14 

63 

56 

70 

57 

79 

18 

85 

24 

68 

61 

72 

58 

26 

91 

31 

72 

65 

74 

58 
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30 

94 

33 

75 

67 

75 

58 

1-00 

.0907 

46 

85 

74 

79 

61 
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1-15 

12 

49 

89 

76 

80 

63 

1-30 

16 

56 

92 

79 

81 

63 
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20 

58 

98 

84 

83 

64 

2-30 

28 

67 

.0702 

86 

85 

65 
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35 

76 

09 

92 

88 

66 

79 

4-19 

44 

84 

14 

96 

91 

68 

4-30 

44 

86 

15 

97 

92 

68 

79 

5-30 

50 

92 

19 

.0500 

94 

68 

6-00 

52 

94 

21 

02 

95 

68 

79 

6-33 

55 

96 

22 

03 

95 

68 
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TEST    21 
4    LB.    PLATE  4.1    psi 


Time 

Temp 

Hrs-Mins 

r=0 

r=4" 

r=l" 

r=l£" 

r=2" 

r=2i" 

°F 

00 

.1610 

79$ 

01 

.1651 

.1572 

.1333 

.0998 

.0377 

.0159 

03 

68 

87 

47 

.1010 

83 

61 

06 

81 

98 

57 

19 

88 

64 

12 

95 

14 

72 

31 

95 

67 

15 

.1700 

19 

77 

34 

97 

69 

79* 

25 

11 

33 

87 

44 

.0402 

70 

30 

14 

37 

90 

47 

04 

71 

80 

51 

27 

48 

.1403 

57 

10 

74 

1-00 

31 

54 

06 

61 

11 

75 

80 

1-30 

41 

63 

16 

68 

16 

79 

2-00 

49 

70 

23 

75 

19 

80 

80 

2-45 

58 

78 

30 

81 

24 

82 

3-00 

61 

81 

32 

83 

25 

83 

80 

3-30 

64 

86 

37 

88 

28 

87 

79$ 

4-03 

69 

88 

41 

91 

31 

89 

4-30 

72 

92 

45 

94 

32 

90 

79- 

TEST 

22 

4  LB. 

PLATE 

4.1  psi 

me 

r=0 

r=£" 

r=l" 

r=l£" 

r=2" 

r=2£" 

Temp 

op 

00 

.1580 

81 

01 

.1610 

.1538 

.1313 

.0973 

.0435 

.0136 

03 

30 

53 

29 

85 

42 

40 

06 

44 

68 

40 

95 

48 

42 

15 

63 

91 

56 

.1011 

56 

46 

24 

75 

.1602 

65 

19 

63 

48 

30 

80 

06 

69 

23 

65 

49 

48 

91 

18 

79 

32 

71 

51 

1-00 

96 

21 

83 

35 

73 

52 

81 

1-30 

.1706 

33 

92 

43 

78 

54 

2-00 

13 

43 

99 

50 

82 

55 

81 

3-05 

25 

54 

.1411 

59 

86 

59 

4-00 

33 

63 

17 

66 

90 

60 

81 

4-30 

37 

68 

20 

69 

95 

62 

5-05 

40 

70 

23 

71 

96 

62 

80i 

6-00 

46 

76 

28 

75 

98 

63 

80 
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TEST 

23 

4  LB. 

PLATE 

6 . 1  psi 

IME 

Temp 

rs-Mins 

r=0 

L     2 

r=l" 

r=l£-H 

r=2" 

r=2^" 

°F 

00 

.1875 

76 

01 

.1918 

.1820 

.1554 

.1138 

.0488 

.0243 

03 

39 

35 

70 

53 

96 

46 

06 

51 

47 

83 

63 

.0503 

51 

15 

66 

66 

95 

77 

12 

54 

25 

75 

78 

.1601 

85 

18 

55 

76 

37 

82 

85 

07 

91 

21 

57 

45 

85 

89 

09 

93 

23 

58 
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89 

94 

13 

94  . 

26 

60 

76 

1-15 

93 

94 

17 

99 

27 

61 
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95 
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18   , 

.1202 

29 

61 

2-00 

99 

02 

23 

05 

30 

62 

76 

3-00 

.2009 

07 

27 

09 

34 

63 

76^ 

4-00 

14 

08 
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12 

35 

64 

77 

5-00 

17 

18 

39 

17 

38 

71 

5-33 

20 

19 

39 

18 

39 

71 

11\ 

TEST    24 


4  LB. 

PLATE 

6.1  psi 

Temp 

me 

r=0 

r=£" 

r=l" 

r=l£" 

r=2" 

r=2fc" 

op 

00 

.1925 

84 

01 
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.1873 

.1583 

.1187 

.0451 

.0212 

03 
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97 

.1605 

.1203 

62 

18 

06 

16 

.1908 

20 

14 

71 

22 

08 

24 
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27 

20 

76 

24 

15 

45 

37 

43 

33 

84 

28 

84 

26 

64 

60 
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47 

95 

32 

30 
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62 

62 
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33 

84 

45 

81 

72 

73 

58 
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90 

81 

81 
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07 

38 
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.2103 

93 

93 

73 

15 

41 

2-00 

14 
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.1702 

80 

20 

44 

84 

2-34 

24 

10 

10 

85 

25 

48 

3-00 

29 

19 

15 

89 

27 

49 

84 

3-30 

35 

20 

20 

94 

30 

50 

4-02 

40 

29 

26 

.1300 

31 

51 

84 

4-30 

45 

33 

29 

03 

34 

52 

5-00 

49 

39 

34 

06 

36 

54 

5-30 

54 

44 

37 

09 

38 

55 

84 
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APPENDIX  G 
Spectrochemical  Analysis  Report 

The  reported  analysis  was  made  by  Mr.  Walter  W.  Correia 
of  the  STRNAD  Spectrographic  Laboratory  at  M.I.T.   The  sample 
labels,  "4",  "6",  and  "8"  are  for  the  4  lb.  plate,  6  lb.  plate, 
and  8  lb.  plate  respectively.   The  analysis  samples  were 
taken  from  tensile  test  specimens  randomly  selected  from 
those  actually  used  in  the  tensile  experiments.   From  this 
data  one  may  infer : 

4  lb.  plate  Pb  >  99.67% 
6  lb.  plate  Pb  >  99.96% 
8  lb.  plate       Pb  >  99.29% 
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MASSACHUSETTS     INSTITUTE    OF     TECHNOLOGY 
STRNAD   SPECTROGRAPH  I C    LABORATORY 
Rm.     13-4151 

Report  of    Spectrochemicol     Anolysis 


To.      Pftt.ftr  T.   Tarpgaard 


Sample  Description: 


Pb 


Charged  To- 

Date   Received^ 

Date   Reported:    March  ?>   ^TO 


SAMPLE 

4          6          8 

SAMPLE  I      4          6          8 

SAMPLE 

4        6            8 

Ag 

VVFT 

VFT 

VFT 

Mg 

VFT1 

VFT 

VFT 

Ti 

Al 

Mn 

Tl 

As 

Mo 

U 

Au 

Na 

V 

B 

Nb 

w 

Ba 

Ni 

Y 

Be 

Os 

Zn 

Bi 

T 

FT-T 

T 

P 

H 

H 

H 

Zr 

C 

Pb 

Ce 

Co 

WFT 

VFT 

WFT 

Pd 

Dy 

Cd 

T 

VFT 

Pt 

Er 

Co 

Rb 

Eu 

Cr 

Re 

Gd 

Cs 

Rh 

Ho 

Cu 

VVFT 

FT 

VFT 

Ru 

La 

Fe 

Sb 

FT-T 

FT-T 

T 

Lu 

Ga 

Sc 

Nd 

Ge 

Se 

Pr 

Hf 

Si 

VFT 

VFT 

WFT 

Sm 

Hg 

Sn 

FT 

L 

Tb 

In 

T 

Sr 

Tm 

Ir 

Ta 

Yb 

~    ' 

K 

Te 

Li 

Th 

?emnrk<;:         Sn  in  8&TnPle   6  approx.   0.5%-Q,J>% 


Vork  by:_ _Book  No  Spec  8   Page  No.    28 

Sported    by-'  ^^0^>7l-( £f<L?^£i  *, 
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Key:   VVFT  <  O.OOOI  % 

VFT  0.0001  %  -0.001  % 
FT  0.001  %-O.OI% 
T  0.01  %   -0.1  % 
L  0.1%  -1.0% 
M    1.0% -10.0% 
H   >  10.0% 


APPENDIX  H 
Computer  Programs 

I .   Program  for  Clamped  Edge  Case 

The  following  listing  is  for  the  program  used  in  the 
numerical  evaluation  of  the  energy  dissipation  rate  integral 
for  the  clamped  edge  circular  plate.   The  program  consists  of 
a  main  part  and  two  subroutines.   Subroutine  SAMS  evaluates 

the  moment  terms,  M  k   +  M-k_,  at  a  radial  position  x,  where 

r  r     6  0  ^ 

x  =  r/a.   Simlarly  Subroutine  JOES  evaluates  the  extensional 
force  terms,  N  e   +  Nfiefl,  at  a  position  x.   The  main  program 
uses  five  point  Gauss  quadrature  to  evaluate  the  total 
integral  over  x.   The  program  is  written  in  FORTRAN  IV. 

The  inputs  to  the  program  are  n,  c' ,  c',w  /h,  and  a,.   The 
input  format  is  F10.5  which  means  each  input  quantity  is 
alloted  10  spaces  and  if  no  decimal  is  inserted  by  the 
programmer  the  decimal  is  assumed  to  be  to  the  left  of  the 
last  five  spaces  in  the  field.   The  format  for  the  input  card 
is  as  follows : 

Field  Quantity 

1-10  n 

11-20  c| 

21-30  c'2 

31-40  w  /h 

o 

41-50  a. 

The  last  data  card  must  have  a  negative  number  in  the  field 
1-10.   This  causes  the  program  to  exit.   For  each  data  card 


118 


the  program  computes  and  prints  seven  sets  of  output  data 
where : 

starting  with  the  input  value  of  c'  plus  .1.   All  other 
parameters  remain  constant  through  all  seven  iterations. 
The  output  is  in  ten  columns  as  follows: 

Column  Quantity 

1  n 

2  c^ 

3  C2 

4  w  /h 

5  a, 

6  A  ,  the  integral  of 

m  * 

M  k   +  Mqk„ 
r  r     0  6 

7  A  ,    the  integral  of 

N  e   +  Nnefl 
r  r     0  6 

8  A  +A   =  I 

m   n    c 

9  $ 
10  <Dn 


where   ^  = 


3-a, 
c 


119 


r\i 
»— • 
Q 

* 

•n 

00 

no 
*o 

CM 
C> 

m 
cm 


UJ 

lO 

< 

H 

u 

< 

Ui 

x 

o 

3: 

o 

» 

UJ 

< 

CM 

Q 

U 

UJ 

« 

Q. 

< 

■57 

.-• 

o 

< 

U 

— 

04 

_1 

• 

LT\ 

* 

O 

2 

• 

O 

LU 

O 

r\l 

^» 

r-l 

* 

O 

U 

O 

vO 

m 

.—« 

UJ 

< 

Q 


rvj 


ii 


O 

• 
* 

UJ 

+ 
< 


X 

> 


UJ 

x 

2 

* 

* 

i-H 
U 


to  X 
5  >: 


x 
> 


UJ 

•> 

I 

"5 

» 

O   .-H 

<f  • 
co  z> 
c> 
r- 


x 

> 


UJ 

I 


u 


—        o  ~ 


rg 


Q 

< 

UJ 

or 

in 


—  cm              <-•  o  <   > 

»-  i  •  u  o  ->  U  •  lO 

<  Uj        cm  oh 

5     -     II     r<     ||      II  _j 

a  ii  _i  c 

OU.IMOTH  <t  *— 
u.— 'UOujuduo 


1/7 
> 
< 


O^   t0 

M     _J 

Z>  W 


O  3 


lO 
< 


X 

> 


UJ 


o  u 


o 

en 


X 

> 


UJ 


osl 


m  to 


ii 


J  <M    H    J  W 

u  o  z>  u  o 


» 

—  z> 

o  ~ 

• 
r-*    lO    X 

—  <  > 

*    CO 

Z>  II 

-J 

II    _J  <f 

<  •— 

Z>  U  O 


UJ  ►-• 
CC  Q 
3  * 
t~  <-* 

<  in 

0£  30 
Q    TO 

<  M3 
3  <M 

O  o 

vO 

</>  en 

tO   (VI 

D    • 

<  -f 

iD  O 

V-  o 

a  * 
a 

m  ao 

30 

r  to 

I-   30 

»-«  ao 

H   CO 

00 

_J    vO 

<  m 

<Y.     • 

<3  — 

UJ 


a 
* 
»n 
o 

r- 

vO 
CO 

f\J 

vO 

CO 

r- 

• 
+ 

CO 


H-  O 
71  m 


UJ    II 


* 
m 
o 
r- 

vO 
CO 
CM 
vO 
CO 

r-  o 

•  o 

ii 


o 
MT> 


O 

CM 


u 


X 

> 


UJ 

x 

CM 


X 

> 


UJ 

X 

cm 


KJ 


ON  t-« 

Z>        OOD 

r- 

uj  2:  vo  uj  r 
o  >  o  o  > 
&>  -) 
•        II 

_» 
II  _J  t-l 
<  — • 


•  -> 


X 

> 


UJ 

X 


C\J 


o  — 

• 

^i  to 

I     UJ 

—  o 
*  -) 


X 

> 


UJ 

» 

I 

CM 


o  u 

ox  — 

<t  to 

ao  UJ 

m  o 


II 


>  in 
II    • 


_J  o 


_l 

_J   <M 


Z)  U  O  3  U  O 


< 


II  -J 
< 
c  z>  u 


120 


CM 

»-« 

a 

* 

•-4 

m 

CO 

CO 

o 
f\J 

:n 

f\| 


O 

h- 

* 

> 

!-• 

* 

in 

— 

f- 

CO 

-t 

> 

CO 

•— ♦ 

m 

vO 

Q 

Z 

— 

CM 

* 

> 

X 

O 

m 

» 

Z 

vO 

o 

5: 

> 

en 

r- 

> 

• 

r\j 

vO 

• 

.-« 

• 

CO 

t-* 

< 

+ 

(M 

< 

• 

O 

vO 

•> 

H 

fr-4 

CO 

--» 

X 

UJ 

O 

r- 

»- 

3: 

•> 

* 

<f 

> 

* 

— 

X 

& 

• 

* 

m 

UA 

>c 

CO 

+ 

o 

U 

• 

•> 

JO 

m 

• 

» 

m 

eg 

CO 

»-■* 

c\j 

r~t 

l-l 

u 

o 

n 

t~« 

V 

Ul 

•» 

CO 

* 

— 

*> 

iTi 

r-l 

co 

m 

\ 

*z 

« 

U 

00 

o 

— 

"T 

Ul 

in 

» 

o 

r- 

r-4 

IU 

— • 

• 

™ 

T) 

<n 

o 

s 

< 

* 

o 

o 

O 

— 

• 

co 

> 

1 

* 

o 

.--» 

• 

— 

f\j 

o 

<f 

r- 

u. 

X 

rH 

IT) 

* 

*o 

+ 

• 

h- 

« 

u\ 

♦- 

z 

1 

UJ 

X 

o 

en 

ro 

>  c 

■— » 

UJ        »-• 

> 

—. 

o 

21 

in 

f- 

> 

—* 

•— 

D  UA    X 

* 

-> 

> 

• 

-* 

> 

II 

y- 

Z         UJ 

N 

3 

li 

• 

M 

UJ 

< 

—  C 

_j 

n 

+ 

n 

■  r 

t— 

5 

»■ —  1 —  1 

en 

II 

_i  ^» 

< 

< 

•— « 

a 

z      _j  a 

»— « 

< 

*-« 

2* 

»- 

l— 

K 

cr 

o 

o  o  <  z 

Q 

D 

u  c 

> 

> 

> 

> 

2E 

u 

U  O  U  UJ 

r-« 

o 
c 
r- 

o      o 

i-t 

121 


* 

■**-». 

O 

• 

CM 

* 

* 

—   CM 

—    < 

O   — 

•    * 

— » 

CM  — 

r 

*    O 

3s  — 

*      • 

*  j: 

— 

—    CM 

<M    jZ 

-^ 

O    * 

U    * 

"NJ 

•    * 

*    -~ 

< 

CM   X 

—    O 

u_ 

*    — 

— »      • 

'1 

Ll   — 

*     * 

O    CM 

1    rg 

—  X  O 

•    * 

CNJ   < 

O    1      • 

ry|    * 

<   * 

•    O  CM 

— ■ 

*     ~ 

UJ  O- 

CM     •   — 

cm 

*     O    -~ 

u.     •  f\i 

*     H     + 

< 

x   •  r 

—  (N  < 

* 

<-» 

* 

—  CM  "2 

*    +  * 

X 

X 

O 

*  *  * 

en  —  0 

*   *  0 

5L 

• 

m  *    — 

lL  —     • 

in  —    • 

> 

0 

•  X  O 

—  —  en 

•   O  CM 

* 

~- 

tn   |     • 

1   0  — 

CM     •   * 

<-* 

+ 

1    O  r\j 

— •    •  -f 

1     CM    * 

< 

— » 

X     •    * 

—  r\j  — 

X    *    X 

•> 

— » 

*      rH     * 

r-4     *      ~~ 

*     *     1 

zr 

— 

m  —  x 

— - 

» 

<     *      — 

O    rH    O 

LU 

O 

•  —    I 

T 

rH 

r\l 

LL    X  3 

— - 

•    <       • 

*> 

• 

-^ 

<M  *    O 

35 

to 

tO 

Ll   —     • 

< 

CM   —    •-* 

r 

CM 

< 

■_  «-.     , 

* 

* 

* 

+    *     CM 

< 

—    *     — 

-5 

4c 

< 

~  O  — « 

— 

r~. 

-~ 

rH    O     * 

LU 

* 

•> 

* 

LL 

+     •  — 

O 

— 

-~ 

<       •     * 

* 

+   O  CM 

<M 

X 

* 

~  CM  * 

• 

— 

;vr 

-» 

i? 

lu  m  x 

< 

~     •  < 

u 

~- 

< 

r  *  cm 

CM 

Z 

LU 

ZT 

Ll 

Ll    1    ~ 

LU 

X    CM    * 

* 

* 

Lu 

3  *  < 

* 

LU 

V. 

LU 

\ 

"  O  * 

* 

jC     *     rH 

-H 

O 

* 

*  -•  * 

* 

V. 

O 

\ 

O 

*     •  O 

— » 

*      *     < 

c; 

• 

-» 

rH    <    rH 

CM 

O 

• 

O 

• 

-~    rH      • 

O 

rH     X      * 

•. 

r- 

O 

u  —  < 

< 

• 

rH 

• 

r-i 

^  ~  m 

• 

\J 

D 

1 

• 

*  *  * 

— ■ 

.— I 

+ 

i-« 

-*- 

<  *    1 

rH 

*   *  0 

w 

• 

rH 

+ 

x  0  r> 

* 

CM 

+ 
O 

O 

• 

<f 

+ 
O 

O 

• 

H.HO 

*  <    • 

+ 

—  r    • 

X    S  CM 

to 

en 

^ 

*    •    • 

O 

rH 

• 

cm 

.-1 

• 

r\l 

0  *  en 

LU 

*   *  * 

V 

1 — 

— » 

LU 

m  (\i  oj 

• 

•> 

rH 

^~ 

^ 

rH 

~- 

•  0  — - 

* 

00* 

< 

* 

< 

* 

•  *    * 

CM 

rvi 

— » 

«-» 

* 

^ 

— 

— 

* 

CM     •    * 

« 

O 

•     •  X 

tO 

.-< 

rH 

'X 

O 

CM    *     * 

* 

r-i 

Ill 

* 

* 

r~\ 

ll 

* 

* 

—     <f     rH 

<t 

• 

(M   <XJ  ■— 

m 

< 

< 

■— 

• 

1     X   X 

* 

< 

* 

U. 

* 

< 

•1 

:l 

* 

< 

\    --    <  LU 

CM 

1     —    * 

< 

'U 

• 

1 

* 

<M 

in  —  — 

vr 

< 

0 

«-H 

«-» 

< 

LU  O 

cn 

--- 

< 

Ll 

O   —   * 

«— 

~^ 

m   +  X 

< 

a 

i' 

+ 

O 

to 

~^ 

•  *    * 

— 

Ll 

U 

.—» 

0 

1/) 

UJ 

Ll  Ll 

r-« 

n 

to 

Ll 

U. 

•    *    O 

1/) 

«— 

•  —  3 

Ll1 

UJ 

•— 4 

D 

O 

• 

n 

•^ 

f-*  O  0 

* 

* 

+ 

0 

0 

• 

m 

Ll 

—     1 

• 

en 

Ll 

*^ 

r-t  ;>     • 

an 

V. 

rH     Z     * 

* 

+ 

(— 

* 

• 

CM 

<t 

>v 

—    •    • 

3 

m 

U 

• 

• 

— 

rH 

< 

U 

— 

•-* 

< 

~-   C£  <M 

<r 

Z 

—   3  O 

m 

U 

D 

tT\ 

<—* 

— 

7* 

—    vO    vO 

• 

• 

U 

.-« 

,--> 

Ui 

1 

II 

II    Ll 

U. 

1 

ll 

II 

^^         w- 

LU 

—  *      . 

• 

LU 

O 

• 

11 

UJ 

c-i  ,-i 

<r 

u 

11 

u 

11 

II 

11 

CM   00   — 

nr 

II 

II 

II 

11   —  — 

— 

11 

11 

II 

II 

— 

II 

rH 

<M    II 

~-~ 

11 

<— « 

CM 

II         II 

II 

11  u-r 

11 

11 

en 

II 

< 

♦  -f 

+ 

< 

< 

< 

< 

< 

< 

V 

< 

*   +  :< 

w-4 

D 

C\i 

< 

< 

a' 

u 

0 

ill 

rH 

rvj 

u. 

rH 

LU 

L11 

LU  Ll 

Ll 

c\. 

u. 

Ll 

LL 

5    5    <- 

*a 

DC 

u 

O 

Ll 

U) 

X 

< 

LL 

LL 

u_ 

Ll 

f~»   CM 

en 

U 

u 

to 

tO 

•— • 

CO 

H 

li 

CM 

U 

Ll  Ll 

en 

u 

Ll 

U 

ck  a:  uj 

UJ 

UJ 

UJ 

rH    CM    CO 

LU 

LU 

LU 

r  1 

*-t 

r— • 

I— 1 

O 

< 

Q 

122 


c\j 

< 

U. 

LU 

1 

f\J 

< 

UJ 

UJ 

w       ^-* 

*   rvj 

03  < 

UJ    * 

■—  O 

1      • 

—  rsi 

—    4- 

H    «-» 

X 

<    — 

* 

T-t 

f\J 

u.  o 

-~ 

00 

I/) 

UJ      • 

— 

* 

* 

+   CM 

~-> 

~ 

-«. 

<~t    $ 

V 

— 

— » 

<    * 

j> 

*— 

?r 

— - 

z 

UJ    X 

f- 

i* 

UJ 

4- 

Ui 

UJ     f 

w 

UJ 

\ 

LU 

\ 

—  a 

;/> 

^ 

o 

-N. 

o 

*    • 

a 

o 

• 

o 

• 

—    r-« 

<r 

• 

r-* 

• 

rH 

— .     »_ 

— ~ 

r-H 

+ 

r-l 

+ 

<     * 

■«- 

+ 

o 

+ 

o 

UJ     r-l 

— 

«c 

o 

• 

o 

• 

*      < 

~~ 

<-l 

• 

CM 

oo 

• 

CNJ 

o  * 

V 

» 

.— « 

— 

<-* 

— » 

.—« 

— * 

•  o 

s: 

vO 

— 

"— 

* 

•k 

u 

w 

* 

r\i     • 

or 

.—1 

UJ 

* 

* 

cd 

UJ 

* 

* 

—   4" 

— 

•• 

UJ 

* 

< 

.— « 

»-. 

* 

< 

\  — 

>T* 

in 

«-» 

< 

UJ 

a 

•» 

<t 

u 

o  — 

11 

r-H 

a 

1/5 

UJ 

u. 

u. 

r- 

o 

LO 

U_ 

Ui 

•  * 

<t 

• 

CO 

LU 

— 

1 

rH 

• 

n; 

til 

~- 

•H    > 

^-> 

o 

q 



r-< 

■a 

u 

— 

rH 

< 

—   I— 

~- 

• 

• 

U! 

1 

11 

II 

u 

Ll 

1 

II 

II 

■■— 

<?* 

r  « 

r  4 

U! 

ii 

LU 

ii 

II 

II 

or 

II 

It 

— 

II 

t-A 

ru 

II 

~- 

II 

.-- 1 

-M 

II 

"> 

<i 

<. 

< 

*-l 

<l 

< 

V 

X 

f^  Q 

J           .-« 

(\J 

It 

•--» 

UI 

UJ 

UJ 

u 

u 

f\J 

u 

u. 

u. 

I-    > 

> 

uj  zr 

\/l 

»/l 

+—• 

f~l 

UJ 

UJ 

UJ 

u 

*v-« 

r-l 

u; 

a 

r     i 

u. 

UJ 

r-    r- 

> 

rr  U! 

t 

123 


x      — 
5      o 


* 

• 

X  — 

rM 

-^ 

<M 

-,  ~ 

"5  ^. 

u 

-~ 

* 

x  — 

*    O 

* 

0 

* 

3  0 

CM     • 

• 

< 

0 

• 
M 

< 
* 

• 

(M 
* 
* 

O 
• 
cm 
* 
* 

X 

• 

+ 
r 

*    • 

CM   CM 

*  * 

O    O 

•      ■ 
CM   (M 

*  * 

CM 
< 

U  CM 

*  * 
—    * 

::>  — 
•  0 

CM     • 

*  cm 

*  * 
X    * 

I! 
< 

1 

X 

~- 

* 

* 

X 

< 

X   X 

0 

* 

—  X 

* 

* 

X 

O 

ro 

1 

* 

*    1 

• 

O 

*  — 

— - 

CM 

z 

• 

— 

O 

— - 

0  0 

CM 

• 

u>   1 

— - 

U 

> 

en 

1 

• 

0 

•      • 

* 

CO 

•  0 

O 

* 

« 

— ■ 

— 

^H 

• 

en  r-i 

* 

^' 

rM     • 

• 

— 

r-i 

+ 

X 

~» 

cm 

1    — 

X 

•1- 

1      r-» 

CM 

— ~ 

< 

— 

* 

— 

* 

X  — 

1 

-^ 

X  — 

* 

O 

•> 

— . 

Lf\ 

* 

* 

*   * 

O 

— » 

*      — 

* 

• 

X 

2: 

— . 

• 

~ 

X 

0  — 

• 

-^ 

^~ 

O     * 

X 

— 

CM 

* 

UJ 

0 

M 

O 

1 

•  0 

.— 1 

X 

u 

r> 

•   X 

»— 

r 

*     U_ 

— 

* 

• 

-^ 

• 

O 

CM     • 

— 

5 

CD 

• 

CM  3 

1 

~2 

*    <J 

~ 

X 

cm 

^- 

CM 

• 

—    CM 

* 

* 

* 

CM 

—  * 

O 

* 

X    * 

— 

3 

* 

— 

* 

.—1 

—    * 

— 

— 

— 

* 

^  .-, 

• 

— 

—  — 

UJ 

» 

* 

+ 

* 

— 

+     * 

O 

O 

Z 

* 

+  O 

r-l 

0 

1    z 

z 

CM 

X 

~ 

r~ I 

* 

. — 

~ 

—    H 

• 

• 

< 

X 

—■     • 



• 

~ 

— 

X   < 

r- 

U 

—* 

i: 

< 

(\J 

z 

z 

X  < 

fM 

rM 

UJ 

— 

x  ™ 

* 

fM 

z 

Z 

—    UJ 

— « 

• 

* 

3 

w 

< 

UJ 

uj 

JS    — ' 

* 

* 

a) 

* 

:-;  * 

CM 

* 

UJ 

UJ 

—  V 

to 

rH 

0 

* 

* 

* 

\ 

\ 

*    * 

* 

* 

1 

O 

*  * 

<1 

* 

V 

\ 

+   1 

CO 

O 

• 

.-1 



r-« 

0 

0 

r-l    ~ 

X 

CM 

2: 

• 

H   i-H 

* 

'M 

0 

0 

—  z 

< 

* 

r-~ 

w 

O 

< 

• 

• 

u  0 

— 

< 

< 

iT» 

u  < 

r^ 

< 

• 

• 

X  < 

— 

D 

1 

* 

• 

* 

r- 1 

.-1 

*         • 

* 

>-~ 

Q 

1 

*  — 

< 

^> 

r-l 

r-l 

3  a 

+ 

— 

0 

— 

'M 

-— 

+ 

+ 

~-   CM 

r 

* 

ca 

O 

—  * 

^ 

* 

+ 

+ 

*  u 

~ 

• 

X 

* 

O 

~ 

0 

0 

X    * 

< 

--^ 

»— 

s: 

• 

X  — 

-~ 

-^ 

0 

O 

r-l    — 

-^ 

to 

en 

* 

* 

• 

1 

• 

• 

*     * 

* 

O 

* 

QZ. 

CO 

*  0 

O 

O 

• 

• 

u  * 

UJ 

«jj 

~~ 

in 

X 

cm 

~^Z 

-^ 

r-l 

r-l 

O  X 

rM 

• 

-» 

UJ 

*— 

0    • 

• 

• 

r-l 

f-» 

*  — 

z 

O 

!— 1 

* 

• 

— • 

* 

* 

a 

— 

-~ 

•«-» 

•    -^ 

< 

CM 

'i 

J— 

* 

•  CM 

CM 

CM 

— 



--- 

«-* 

—  < 

Of 

~) 

< 

rH 

fM 

* 

* 

— « 

no 

* 

UJ 

A 

CM   * 

* 

* 

GC 

r-l 

CM   * 

* 

* 

O 

* 

UJ 

* 

x  u 

— 

if\ 

< 

1 

I 

X 

0 

^~ 

* 

CO 

* 

1     X 

.— j 

* 

\ 

'3C 

< 

1     * 

* 

* 

u 

* 

*^» 

* 

1     V; 

lO 

UJ 

• 

1 

* 

m 

^ 

— 

• 

< 

W 

< 

u 

— 

< 

O   2: 

< 

X 

O 

liJ 

* 

m  x 

X 

X 

< 

u 

" 

< 

w 

— 

< 

O    O 

CO 

z 

+ 

O 

• 

* 

* 

CM 

CC 

QC 

m 

a 

QC 

1/1 

Uu 

•   * 

* 

— 

m 

* 

O 

•  — 

«— 

~- 

u 

u 

I/O 

O 

u 

^ 

UJ 

•      • 

< 

>— 1 

J 

O 

• 

t-H 

O 

O 

* 

* 

+ 

ao 

CO 

1 

GO 

03 

t-H  O 

O 

* 

i-H 

or 

• 

r-«    * 

* 

* 

* 

+ 

<x> 

KJ 

1 

CD 

u 

r-l    r-l 

1— 

+ 

• 

CM 

— 

• 

• 

* 

a> 

ce 

< 

(T 

< 

~-»     • 

• 

O 

— 

z 

CM 

*-    O 

D 

CD 

m 

X 

< 

a~ 

< 

-^    

D 

ir 

."1 

— 

— 

^O 

nO 

C\, 

• 

on 

11 

or 

11 

~    vf> 

vO 

• 

— 

—       • 

• 

• 

• 

u 

II 

U 

II 

~- 

O 

• 

<~i 

r-« 

<L 

11 

11 

1— 1 

«— 1 

<r 

II 

10 

00 

a> 

CM 

II 

II 

II 

11  0 

C* 

It 

II 

II 

~~. 

— 

■w 

11 

11 

z 

11 

Z 

II   — 

— 

-^ 

>— « 

II 

II     — 

^-r 

— 

11 

II 

Z 

11 

z 

II 

Ct 

11 

+ 

f 

* 

<t 

< 

< 

<l 

+ 

+ 

+ 

UJ 

+ 

+ 

+ 

< 

< 

< 

< 

UJ 

r-l 

13 

f\i 

< 

a 

KJ 

0 

a 

0 

U! 

UJ 

UJ 

u. 

ir. 

m 

< 

CC. 

u 

Q 

0 

a 

Ui 

Uj 

UJ 

u.  z 

ui 

X 

<T 

<r 

a" 

(I 

CO 

a 

cc: 

a' 

a 

or 

cr 

a 

*— 1 

U 

u 

w 

U 

(J 

0 

u 

^ 

u 

O    r- 

•~t 

cm 

CO 

•—• 

cv 

CO 

2 

1— 

r-l 

CM 

CO 

UJ 

<£ 

< 

Q 

u 

124 


II.   Program  for  Simply  Supported  Case 

The  following  program  for  the  simply  supported  case  is 
exactly  like  that  for  the  clamped  case  except  that  a  different 
assumed  velocity  profile  is  used. 

The  input  data,  in  format  F10.5,  is  as  follows: 
Field  Quantity 

1-10  n 

11-20  c' 

21-30  c'2 

31-40  w  /h 

Here  again  the  last  data  card  must  have  a  negative  number  in 
field  1-10  to  cause  the  program  to  exit.   The  output  data  is 
exactly  as  in  the  clamped  edge  program  discussed  in  Part  A 
except  that  the  column  for  a,  is  omitted  and  here 
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APPENDIX  I 
Computation  of  Total  Deflection 

The  discussion  in  Chapter  II  was  centered  about  determining 

the  deformation  rate,  w  .   In  practical  problems  it  is 

frequently  the  total  deflection,  w  ,  at  some  time,  t* ,  which  is 

w 

of  interest.   It  has  been  shown  that  w  is  a  function  of  r— , 

o  h 

and  also  that  the  material  creep  constants,  K  and  n,  may 

change  with  time.   To  compute  w  ,  therefore,  a  piecewise 

linearization  technique  can  be  employed  in  which  the  deformation 

process  is  broken  into  a  sufficient  number  of  linearized  time 

segments  to  provide  the  accuracy  desired.   One  such  technique 

is  discussed  in  this  Appendix. 

In  dividing  the • deformation  process  into  time  segments 

some  guidance  can  be  obtained  from  the  shape  of  the  tensile 

w 
test  curVes  and  from  the  shape  of  the  $  versus  r- —  curves , 

Figures  2.2  and  2.4.   The  time  segments  should  be  relatively 

short  during  primary  creep  and  may  be  longer  when  deformation 

is  proceeding  more  slowly.   As  the  deformation  process 

w  w 

proceeds,  r— -  increases;  and  when  r— -  has  changed  enough  to  make 

a  significant  change  in  $  a  new  linear  segment  should  be 

introduced. 

Having  divided  the  total  time  of  interest,  from  t=0  to 

t=t*,  into  n  segments  one  can  proceed  as  follows: 

Let  t=t   be  the  time  at  the  beginning  of  the  segment 
under  consideration  and  t=tf  the  time  at  the  end.   From 
tensile  test  data  determine  the  appropriate  values  for 
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the  creep  constants,  K  and  n,  at  the  time  t  ,  in  the 
center  of  the  given  time  segment,  where: 

t  -t 

t   =  -£=-£.  (I.D 

c     2 


Using  these  values  of  K  and  n;  and  using  the  value  of 

w 

7—  at  the  beginning  of  the  time  segment,  compute  w 

with  the  appropriate  equation  from  Chapter  II,  (2.53)  for 

clamped  edge  or  membrane  or  (2.57)  for  simple  support. 

Multiply  this  value  of  w  by  the  total  time  length  of 

the  segment,  tf-t  ,  and  add  the  result  to  the  value  of 

w   at  t  .   This  gives  w   at  t,.   The  procedure  is  then 

os        3      o     f 

repeated  for  the  next  segment,  and  the  next,  and  so  forth 

until  t*  is  reached. 

In  illustration,  consider  the  case  of  Test  13,  a  6  lb. 
plate  under  6.1  psi  pressure.   Assume  that  it  is  desired  to 
compute  w   at  t=8  hours.   The  total  time  is  broken  into  two 
segments,  a  relatively  short  initial  segment  from  t=0  to  t=2 
hours,  and  a  longer  segment,  where  the  deformation  is  more 
gradual,  from  t=2  hours  to  t=t*=8  hours.   Values  of  K  and  n 
are  determined  at  the  median  times  of  the  segments,  which  in 
this  case  are  t=l  hours  and  t=5  hours.   From  the  tensile  test 
results  in  Chapter  II  the  appropriate  values  are  seen  to  be: 

t=l  hour       n=5.0        K=4.1  x  10~18  (1.2) 

t=5  hours      n=5.0        K=1.7  x  10~18  (1.3) 

In  test  13  the  initial  deflection  was  .0880"  giving: 

w 

r^   =  .840.  (1.4) 

h 
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If  test  data  is  not  available  the  initial  deformation  can  be 

computed  using  elastic  or  plastic  analysis.   Using  (1.2)  and 

(1.4)  and  the  appropriate  value  of  $  from  Figure  2.2  in 

equation  (2.48)  the  value  of  w   for  the  first  segment  is  found 

to  be  .0090"/hr.   Deflection  at  t=2  hours  may  now  be  determined: 

(wq)     =  .0880  +  2(.0090)  =  .1060  (1.5) 

which  gives: 

w 
(~)t=2  =  1.01  (1.6) 

Using  (1.3)  and  (1.6)  with  Figure  2.2  and  equation  (2.48)  we 
find  for  the  second  segment  that  w  =  .0020"/hr.  Deflection 
at  t=8  hours  is: 

(w  ).  0  =  .1060  +  6(.0020)  =  .1180"  (1.7) 

O  t —  o 

A  plot  of  this  two  step  linearization  and  of  the  measured 
experimental  curve  is  shown  in  Figure  1.1.   Also  shown  is  a 
three  step  linear  approximation  made  with  the  technique 
described  above  with  periods  0  <  t  <  1  hour,  1  <_  t  _<  3  hours, 
and  3  <_   t  <_  8  hours.   Still  closer  approximations  to  the 
actual  curve  could  be  made  by  dividing  the  time  period  into 
more  segments . 

In  the  practical  application  of  these  results  it  is 
important  that  initial  deflection  of  the  plate  be  known. 
Computing  this  initial  deflection,  however,  may  be  no  simple 
matter.   Finite  deflection  analysis  techniques  for  plastic  and 
elasto-plastic  deformation  of  plates  are  usually  quite  complex. 
For  materials  in  which  creep  is  important  there  is  often  a 
high  degree  of  sensitivity  to  the  speed  of  application  of  load 
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[28]  which  makes  even  the  selection  of  appropriate  elastic 
and  plastic  material  constants  a  rather  difficult  task.   In 
fact,  such  elasticity  concepts  as  yield  stress  and  Young's 
modulus  may  become  rather  meaningless  for  some  materials  under 
conditions  where  creep  is  important. 

In  cases  where  the  initial  deflection  is  expected  to  be 
quite  small  and  elastic,  one  could  calculate  initial  deflection 
using  the  bending  solution  of  Timoshenko  and  Woinowsky- 
Krieger  [29].   For  larger  elastic  deformation  the  large 
deflection  solution  of  Timoshenko  and  Woinowsky-Krieger  [30] 
is  available.   The  large  deflection  elasto-plastic  deformation 
of  circular  plates  is  treated  by  Ohashi  and  Murakami  [31] [32] , 
and  the  large  deflection  of  elastic-perf ectly  plastic  circular 
plates  has  been  analyzed  by  Crose  and  Ang  [33] . 
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